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ABSTRACT 

We present two-dimensional hydrodynamic simulations of stellar core collapse and develop the frame- 
work for a detailed analysis of the energetic aspects of neutrino-powered supernova explosions. Our 
results confirm that the neutrino-heating mechanism remains a viable explanation of the explosion 
of a wider mass range of supernova progenitors with iron cores, but the explosion sets in later and 
develops differently than thought so far. The calculations were performed with an energy-dependent 
treatment of the neutrino transport based on the "ray-by-ray plus" approximation, in which the neu- 
trino number, energy, and momentum equations are closed with a variable Eddington factor obtained 
by iteratively solving a model Boltzmann equation. We focus here on the evolution of a 15 Mq pro- 
genitor and provide evidence that shock revival and an explosion are initiated at about 600 ms after 
core bounce, powered by neutrino energy deposition. This is significantly later than previously found 
for an 11.2 Mq star, for which we also present a continuation of the explosion model published by 
Buras et al. The onset of the blast is fostered in both cases by the standing accretion shock instability 
(SASI). This instability exhibits highest growth rates for the dipole and quadrupole modes, which 
lead to large-amplitude bipolar shock oscillations and push the shock to larger radii, thus increasing 
the time accreted matter is exposed to neutrino heating in the gain layer. As a consequence, also con- 
vective overturn behind the shock is strengthened, which otherwise is suppressed or damped because 
of the small shock stagnation radius. When the explosion sets in, the shock reveals a pronounced 
global deformation with a dominant dipolar component. In both the 11.2 M0 and 15 A/© explosions 
long-lasting equatorial downflows supply the gain layer with fresh gas, of which a sizable fraction 
is heated by neutrinos and leads to the build-up of the explosion energy of the ejecta over possibly 
hundreds of milliseconds. A "soft" nuclear equation of state that causes a rapid contraction and a 
smaller radius of the forming neutron star and thus a fast release of gravitational binding energy, 
seems to be more favorable for the development of an explosion. Rotation has the opposite effect 
because in the long run it leads to a more extended and cooler neutron star and thus lower neutrino 
luminosities and mean energies and overall less neutrino heating. Neutron star g-mode oscillations, 
although we sec their presence, and the acoustic mechanism play no important role in our simulations. 
While numerical tests show that our code is well able to follow also large-amplitude core g-modes if 
they are instigated, the amplitude of such oscillations remains small in our supernova runs and the 
acoustic energy flux injected by the ringing neutron star is minuscule compared to the neutrino energy 
deposition. 

Subject headings: supernovae — hydrodynamics — neutrinos 



1. INTRODUCTION 

The physical processes that start the explosion of mas- 
sive stars are still not well understood, although this is 
of crucial importance for predicting supernova and rem- 
nant properties, nucleosynthesis conditions and yields, 
and the observable signals from supernovae like neutri- 
nos and gravitational waves. Neutrinos are thought to be 
the main agent of energy transport and energy loss from 
collapsing stars, an expectation that was spectacularly 
verified by the detection of two dozen neutrinos in con- 
nection with Supernova 1987A. Colgate & White (1966) 
suggested already 40 years ago that the intense fiux of 
neutrinos radiated from the nascent neutron star might 
deposit the energy needed to reverse the stellar collapse. 
The detailed physics of this neutrino-driven mechanism, 
however, was worked out only later, based on numerical 
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results of Wilson (1985). Bethe & Wilson (1985) showed 
that neutrinos deposit their energy behind the stalled 
supernova shock front mainly by the absorption of v^. 
and Pe on free neutrons and protons, respectively, which 
are the most abundant nuclear species in the postshock 
layer. Provided the neutrino luminosities are sufficiently 
large, this energy input was found to be so strong that 
the stagnating shock can be "revived" and thus accel- 
erates outward to propagate through the overlying, still 
collapsing layers of the star. 

The viability of this neutrino-heating mechanism has 
recently been demonstrated for stars near the lower mass 
end of supernova progenitors for stiff as well as soft nu- 
clear equations of state (Kitaura, Janka, & Hillebrandt 
2006; Janka, Marek, & Kitaura 2007; Janka et al. 2008). 
The investigated progenitor with a main-sequence mass 
of 8.8 Mq (Nomoto 1984, 1987) can be considered as rep- 
resentative of the ^ 8-10 Mq range. The cores of such 
stars consist of oxygen, neon, and magnesium instead of 
iron and possess an extremely steep density gradient at 
their surface. The latter fact leads to a very rapid de- 
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crease of the mass accretion rate onto the forming neu- 
tron star, cnabhng the stalled shock to continue its ex- 
pansion. This creates favorable conditions for neutrino 
energy deposition, thus allowing a neutrino- driven bary- 
onic wind to be launched, which is sufficiently powerful 
to eject the gravitationally loosely bound outer layers 
of the star. The new results confirm qualitatively pre- 
vious simulations by Mayle & Wilson (1988), although 
the more sophisticated treatment of the neutrino trans- 
port in the new models, which were computed with the 
Prometheus- Vertex code, leads to the predictions of 
a lower explosion energy and a lower neutron excess in 
the ejecta. 

In case of more massive progenitors, Wilson & Mayle 
(1988, 1993) could obtain explosions only by assuming 
that the neutrino luminosities and thus the neutrino en- 
ergy deposition behind the stalled shock were boosted 
by enhanced neutrino transport in the neutron star. Ac- 
cording to Wilson & Mayle this could happen, e.g., be- 
cause of the presence of neutron-finger instabilities below 
the neutrinosphere. The existence of neutron-finger un- 
stable conditions, however, was repudiated on groimds of 
a detailed analysis of the neutrino transport conditions 
in supernova cores (Bruenn & Dineva 1996). Neverthe- 
less, numerical simulations (Buras et al. 2006b; see also 
Dessart et al. 2006) showed that Ledoux convection oc- 
curs inside the proto-neutron star and later than about 
200 ms after bounce leads to accelerated lepton number 
and energy losses because of significantly increased muon 
and tau neutrino luminosities. The fluxes of and z>e, 
which are mostly responsible for the shock heating, how- 
ever, hardly change on a post-boimce timescale of some 
hundred milliseconds, and the mean energies of the ra- 
diated neutrinos and antineutrinos of all flavors even de- 
crease compared to models that ignore convection. 

A potentially stronger multi-dimensional effect that is 
clearly supportive for delayed shock revival, was discov- 
ered to be convective overturn in the layer behind the 
stagnant shock. In this region neutrino heating tends 
to create a negative entropy gradient, because the en- 
ergy deposition is maximal just outside of the gain ra- 
dius (Herant et al. 1994; Burrows, Hayes, & Fryxell 1995, 
Janka & Miillcr 1996). Cooler matter that is accreted by 
the shock is carried in narrow downdrafts closer to the 
gain radius, where it readily absorbs energy from the 
neutrinos. At the same time, heated matter with its 
higher entropy becomes buoyant and starts rising and 
expanding, thus reducing energy loss by the rcemission 
of neutrinos and pushing the shock farther out. Con- 
vective overturn therefore improves the efficiency of the 
neutrino energy transfer because it allows more matter to 
be exposed to strong heating near the gain radius. Effec- 
tively this means that the mass in the gain layer grows 
and that gas stays there for a longer time. This is re- 
flected by an increase of the mean timescale Tadv needed 
by accreted matter to move inward from the shock to 
the gain radius. A larger value of Tadv implies that the 
ratio of the advection timescale to the neutrino-heating 
timescale, Thcat, gets closer to unity, which is a necessary 
condition for a model to approach an explosion (see the 
analysis in Buras et al. 2006b; Thompson, Quataert, & 
Burrows 2005; Murphy & Burrows 2008). 

Convection in the flow between shock and gain layer, 
however, is subject to strong damping because the gas 



is falling rapidly and the time for seed perturbations 
to grow is therefore short mitil the gas crosses the 
gain radius and reaches the Ledoux-stable cooling region 
(Foglizzo, Scheck, & Janka 2006). Convective activity is 
therefore not able to develop in the gain layer despite the 
negative entropy gradient there, unless the initial density 
perturbations arc sufficiently large or the ratio of ad- 
vection timescale to convective growth timescale, Tconvi 
exceeds a critical value (for more details, see Scheck et 
al. 2008). Two-dimensional (2D) supernova simulations 
(Buras et al. 2006a,b) with the Prometheus- Vertex 
code, employing a sophisticated energy-dependent de- 
scription of the neutrino transport, indeed show that con- 
vective overturn in the hot-bubble region becomes much 
less vigorous than in previous multi-dimensional models 
with a grey neutrino diffusion scheme (e.g., Herant et al. 
1994: Burrows et al. 1995: Fryer 1999; Fryer & Warren 
2002, 2004), mainly because the accretion shock in the 
new models lingers at rather small radii and the corre- 
spondingly high infall velocities in the postshock region 
are unfavorable for the growth of convection. Scheck et 
al. (2008), performing parametric hydrodynamic studies 
of the post-bounce accretion phase in collapsing stellar 
cores, found that convection is suppressed when (i) the 
neutrino heating is too weak to produce a steep entropy 
gradient behind the shock, or (ii) the rapid contraction 
of the nascent neutron star leads to shock retraction and 
thus causes very short infall timescales of the gas between 
shock and gain radius. 

In such a situation, however, the standing accretion 
shock instability (SASI; Blondin, Mezzcappa, & De- 
Marino 2003), which is a generic hydrodynamic insta- 
bility of the shocked accretion flow to non-radial de- 
formation modes, can exhibit particularly large growth 
rates (cf. Scheck et al. 2008) and becomes the dominant 
multi-dimensional phenomenon to initiate shock expan- 
sion. This was seen in the studies by Scheck et al. (2008) 
and is confirmed by the linear stability analysis of Ya- 
masaki & Yamada (2007). The SASI stimulates large- 
amplitude bipolar shock oscillations, which create strong 
entropy variations in the postshock flow and thus trig- 
ger violent secondary convection. This helps pushing the 
stalled shock to larger radii and therefore stretches the 
time accreted matter is exposed to neutrino heating in 
the gain layer, establishing healthy conditions for the suc- 
cess of the neutrino-driven mechanism. 

Such a decisive role of the SASI was indeed observed 
in case of an 11.2 A'/q progenitor, whose collapse and ex- 
plosion were simulated with the PROMETHEUS- Vertex 
code (Buras et al. 2006b). Only when the setup of the 
numerical grid did not prohibit the growth of the low- 
est (dipolar and quadrupolar) SASI modes^, which have 
the largest growth rates, could an explosion be obtained. 
Convection alone was too weak to yield sufficient support 
for the neutrino-heating mechanism. Instead it was the 
SASI in the first place which pushed the shock farther 
out and helped to bring the ratio of advection timescale 
to neutrino-heating timescale closer to the critical value 
of unity, which is necessary for reversing accretion to ex- 

In previous simulations with failed explosions (Buras et al. 
2003a) the computational volume was constrained to a lateral 
wedge of 90 degrees (±45°) around the equatorial plane, using pe- 
riodic boundary conditions. This setup prevented the occurrence 
of the dipole and quadrupole SASI modes. 
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plosion. 

The underlying physical mechanism that is responsible 
for the SASI phenomenon is still controversial and cur- 
rently a matter of vivid debate. Foglizzo et al. (2007) 
explain the SASI by an advective-acoustic cycle (see also 
Foglizzo 2001, 2002, 2008). This interpretation is in 
agreement with many properties of the instability ob- 
served in the hydrodynamic simulations of Scheck et al. 
(2008) and of Ohnishi, Kotake, & Yamada (2006). In 
contrast, Blondin & Mezzacappa (2006) and Blondin & 
Shaw (2007) advocate a purely acoustic cycle as the cause 
of the low-mode SASI and see this hypothesis supported 
by their numerical studies. 

In this paper we will present evidence that the SASI- 
aided neutrino-heating mechanism may initiate an ex- 
plosion not only in the case of an 11.2 progenitor 
as recently seen by Buras et al. (2006), but in a wider 
mass range of progenitor stars. We study here a 15 Mq 
model (sl5s7b2 from Woosley & Weaver 1995) which 
is frequently used in previous and present core-collapse 
studies because it can be considered as representative of 
stars in a larger mass interval. Stellar evolution calcula- 
tions by Woosley, Heger, & Weaver (2002), for example, 
produced a very similar core structure also for other su- 
pernova progenitors up to about 20 Mq (see model s20.0 
in Figs. A.2 and A.3 of Buras et al. 2006b). We discuss 
results of a 2D simiilation in which a modest amount of 
rotation was assumed and which develops an explosion 
at 600 ms after core bounce, and compare it with non- 
exploding spherically symmetric (ID) simulations on the 
one hand, and 2D models without rotation on the other, 
in which two different equations of state (EoSs) for su- 
pernova matter are used, namely the Lattimer & Swesty 
(1991) EoS with a soft nuclear phase and the Hillebrandt, 
Wolff, & Nomoto (1984; see also Hillebrandt & Wolff 
1985) EoS with a significantly stiffer nuclear phase (Hille- 
brandt 1994). Both are consistent with current neutron 
star mass determinations but the former leads to signif- 
icantly more compact proto-neutron stars than the lat- 
ter. The rotating pre-collapse iron core has a spin pe- 
riod of about 12 seconds. This is roughly a factor of 
ten more rapid than predicted by recent stellar evolu- 
tion models including magnetic fields (Heger, Woosley, 
& Spruit 2005), but considerably slower than necessary 
for magnetohydrodynamics to be important in triggering 
the explosion, which requires pre-collapse spin periods 
of 2 3 seconds (Burrows et al. 2007a, Thompson et al. 
2005). 

We have also c;vahiated our simulations for the pres- 
ence of gravity waves in the nascent neutron star. Bur- 
rows et al. (2006, 2007b) found in their 2D supernova 
models that the anisotropic accretion flow (as a conse- 
quence of the SASI) excited the neutron star to vigorous 
g-mode oscillations, in particular also the dipole [l — 1) 
mode. Later than about one second after bounce, the 
amplitude of these oscillations became very large (about 
3 km) at the neutron star surface and the neutron star 
as a whole showed a large-amplitude periodic displace- 
ment along the z-axis of the numerical grid (Burrows et 
al. 2007c). As a consequence, a sizable acoustic energy 
flux was driven into the surrounding gas and was trans- 
porting a significant amount of acoustic energy to the 
stalled supernova shock. Burrows et al. did not find any 
neutrino-driven explosions in their models. However, at 



late times (typically at t > 1 s after bounce), and there- 
fore much later than the shock is revived by the neutrino- 
heating mechanism in the 11.2 Mq simulation of Buras 
et al. (2006b) and in the 15 model discussed here, 
the acoustic energy flux caused by the core g-mode oscil- 
lations became dominant compared to neutrino heating 
and was able to cause an explosion. In our simulations 
the g-mode oscillations of the newly formed neutron star 
have always such small amplitudes that the acoustic en- 
ergy input to the developing blast is negligible compared 
to neutrino energy deposition. In order to figure out 
whether our neutrino-hydrodynamics code is able to deal 
with this potentially important physical effect, we have 
also performed numerical tests. These demonstrate that 
our code is well able to track large-amplitude g-modes, 
also of dipole character, if such modes are excited in the 
neutron star core. We do not find any reason why our 
third-order hydrodynamics scheme with the chosen nu- 
merical resolution should not be able to follow the long- 
term excitation of neutron star oscillations by anisotropic 
accretion and turbulence in the SASI region between neu- 
tron star and shock front. Whether the excitation could 
be more efficient at times later than currently followed 
in our simulations (at most about 630 ms after bounce) 
remains unanswered and must be left for future explo- 
ration. 

The reason why Burrows et al. (2006, 2007b) did not 
see any neutrino-driven explosions, which is in contrast 
to our findings, is not clear. However, the numerical 
approaches of both groups and the included physics are 
different in many aspects, and therefore it is not aston- 
ishing that for example the rate of neutrino energy de- 
position behind the supernova shock differs significantly 
between the calculations. While our ■■ray-by-ray plus" 
treatment of neutrino transport is highly sophisticated in 
dealing with the energy dependence of the problem but 
describes its dependence on the polar angle only approx- 
imately. Burrows et al. applied a 2D neutrino diffusion 
scheme, in which the energy dependence was only incom- 
pletely taken into account, because energy-bin coupling 
was ignored and Doppler effects due to the motion of 
the stellar fluid were neglected. In addition, the Tucson 
collaboration made use of a different numerical method 
to solve the hydrodynamics equations, e.g. the gravita- 
tional force was implemented in the momentum equation 
in an uncommon, momentum-conserving way. In their 
Vulcan code they employed a polar grid only in the 
outer regions of the computational volume but cartesian 
coordinates in the central part. They performed their 
simulations with the nuclear EoS of Shen et al. (1998), 
which is significantly stiffer than the Lattimer & Swesty 
(1991) EoS, though fairly similar to the Hillebrandt & 
Wolff EoS (Marek et al., in preparation; see also Janka 
et al. 2007). Moreover, their models were computed 
with purely Newtonian gravity, whereas a significantly 
stronger relativistic gravitational potential according to 
Marek et al. (2006) is applied in our models. 

In fact, Newtonian 2D simulations in which the conflu- 
ence of neutrino heating, the SASI, and nuclear burning 
produced explosions of 11 and 15 progenitors, have 
recently been reported by Bruenn et al. (2006) and Mez- 
zacappa et al. (2007) and seem to support the successful 
shock revival found in an 11.2 star by Buras et al. 
(2006b). All these results were obtained with the EoS of 
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Lattimcr & Swcsty (1991). Like Buras et al. (2006b), the 
Oak Ridge group employed for the neutrino transport in 
two dimensions the "ray-by-ray-plus" approximation to 
deal with the lateral dimension. Their transport scheme, 
however, was based on a flux-limited difi^usion descrip- 
tion, in which the energy-dependence of the transport 
was treated in full detail. We consider the discrepant 
results obtained by the Tucson and Oak Ridge collabo- 
rations as additional motivation for us to contribute to 
the discussion by presenting our current simulations in 
this publication. 

The paper is structured as follows: In Sect. 2.1 we give 
a brief summary of the numerical scheme and physics 
input used in the presented simulations, an overview of 
which follows in Sect. 2.2. Section 3 contains our results. 
For reference, we present in Sect. 3.1 results of our ID 
simiilations for the two employed nuclear EoSs, then de- 
scribe the dynamical evolution of our exploding model 
in Sect. 3.2, compare the different computational runs in 
Sect. 3.5, and discuss our neutrino results in Sect. 3.6. 
An analysis of neutron star g-modes including numerical 
tests is presented in Sect. 4. In Sect. 5 we summarize our 
findings and draw conclusions. 

2. METHODS AND MODELS 

In this section we briefly summarize the basic prop- 
erties of our numerical code and the input physics used 
in this work. We then give a compilation of the core- 
collapse simulations whose results are described after- 
wards. 

2.1. Numerical approach and input physics 

The core-collapse and post-bounce calculations 
presented in this work are performed with the 
Prometheus- Vertex neutrino-hydrodynamics code, 
details of which were published by Rampp & Janka 
(2002) and Buras ct al (2006a). The code module that 
integrates the nonrelativistic hydrodynamics equations 
is a conservative and explicit Eulerian implementation 
of a Godunov-type scheme with higher-order spatial and 
temporal accuracy. The self-gravity of the stellar gas 
is treated in an approximation to general relativity, in 
which the monopole term of the Newtonian potential 
is replaced by an effective relativistic scalar potential 
that was constructed from close comparison of the 
Newtonian and relativistic equations of motion (for 
details, see Marek et al. 2006). For this effective 
relativistic monopole potential we consider the Cases A 
and R of the Marek et al. paper. The higher-order terms 
of the gravitational potential are taken as solutions 
of the Poisson equation for the two-dimensional gas 
distribution (see Miiller & Steinmetz 1995). 

The time-implicit transport routine solves the moment 
equations of neutrino number, energy, and momentum 
in spherical symmetry, employing a variable Eddington 
factor for closing the set of equations (Rampp & Janka 
2002). The closure factor is obtained from a simplified 
Boltzmann equation, whose integro-differential character 
is tackled by iterating the coupled systcnn of Boltzmann 
and moment equations until convergence is achieved. 
Our treatment of the neutrino transport and of neutrino- 
matter interactions is energy dependent and includes the 
full redistribution of neutrinos in energy and momentum 
space by scattering reactions with electrons, positrons. 



and nucleons. Effects caused by the motion of the stellar 
fluid are taken into account as well as general relativistic 
redshifting and time dilation. In two spatial dimensions 
the transport is based on the "ray-by-ray plus" approx- 
imation, in which spherically symmetric transport prob- 
lems are solved for the conditions present in the radial 
"rays" corresponding to the different lateral bins of the 
polar grid. This means that we assume the local neu- 
trino phase-space distribution to be axially symmetric 
around the radial direction and thus to depend only on 
one instead of two angles, as a consequence of which the 
neutrino pressure tensor remains diagonal and the lateral 
component of the neutrino flux vanishes. This simplifles 
the transport from a five dimensional, time-dependent 
problem (with variables being the radius, energy, polar 
angle, and two angles that characterize the neutrino mo- 
mentum direction) to only four dimensions and reduces 
the complexity of the moment equations significantly (see 
Appendix B in Buras et al. 2006a) at the expense of neu- 
trino shear and non-radial flux of the neutrinos, which 
are disregarded in such an approach. However, the mo- 
ment equations still contain the terms that account for 
the advection of trapped neutrinos with the stellar fluid 
motion in the lateral direction. Moreover, the lateral 
component of the neutrino pressure gradient is included 
in the hydrodynamics equations. 

Comparison of our approximative treatment of general 
relativity with fully relativistic simulations revealed ex- 
cellent agreement in spherical symmetry and still good 
agreement in two dimensions when the rotation did not 
become too extreme (see Marek et al. 2006, Liebendorfer 
et al. 2005). 

A state-of-the-art treatment of the interactions of neu- 
trinos {v) and antineutrinos (P) of all flavors (Buras et al 
2006a, Marek et al. 2005) is applied in our simulations, 
including all relevant /3-processes, thermal processes, and 
the scattering off electrons, positrons, nucleons, and nu- 
clei. In neutral-current and charged-current neutrino- 
nucleon reactions the effects of nucleon recoil and ther- 
mal motions as well as nucleon correlations in the dense 
medium are taken into account (see Buras et al. 2006a, 
and references therein). Neutrino- nucleon scattering just 
like neutrino-electron scattering is therefore a channel 
for energy exchange between neutrinos and the stellar 
medium and fosters neutrino thermalization, thus lead- 
ing to a noticeable reduction of the mean energy and a 
spectral pinching of radiated neutrinos, in particiflar of 
muon and tau flavor (Keil, Raffelt, & Janka 2003). Elec- 
tron captures on heavy nuclei are described by making 
use of the improved data of Langanke et al. (2003), who 
constructed a table for the capture rates of a large en- 
semble of nuclei in nuclear statistical equilibrium. Also 
interactions of neutrinos of different flavors, scatterings 
as well as neutrino- antineutrino pair conversion (Buras 
et al. 2003b) are included in our treatment. 

We employ two diflferent nuclear equations of state 
(EoSs) in our simulations. On the one hand, we use 
a soft version of the compressible liquid drop EoS of 
Lattimer & Swesty (1991; see also Lattimer et al. 1985; 
"LS-EoS" ) with an incompressibility modulus of bulk nu- 
clear matter of 180 MeV and a symmetry energy param- 
eter of 29.3 MeV. On the other hand we perform simu- 
lations with the considerably stiffer EoS of Hillebrandt 
& Wolff (1985, see also Hillebrandt et al. 1984; "HW- 
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EoS" ) , which is based on a Hartrcc-Fock calculation, as- 
suming a Skyrmc force for the nucleon-nucleon interac- 
tion with parameters given by Kohler (1975). Its incom- 
pressibility has a value of 263 MeV and the symmetry 
energy was chosen to be 32.9 MeV. Details of the cal- 
culation can be found in Hillebrandt et al. (1984) and 
HiUebrandt & Wolff (1985). Both EoSs yield different 
maximum masses of stable, nonrotating neutron stars. 
In case of the LS-EoS the maximum gravitational mass 
is 1.84 Mq, whereas it is 2.21 M© for the HW-EoS, which 
are both compatible with measured neutron star masses 
(Lattimer & Prakash 2007). Both EoSs also lead to a sig- 
nificantly different evolution of the radius of the nascent 
neutron star as a function time (see Fig. 7 in Janka et 
al. 2007). While the LS-EoS produces rather compact 
neutron stars with a radius of i?„s = 11.9 km for a cold 
star with 1.4 and i?ns = 10.0 km for a cold star near 
the maximum mass, the corresponding radii are 13.94 km 
and 13.5 km, respectively, for the HW-EoS (cf. Marek 
2007). 

We apply the nuclear EoS (either LS or HW) only 
above a transition density pEoS- In the low-density 
regime, p < pEoS, the EoS includes the ideal gas con- 
tributions of electrons and positrons (of arbitrary de- 
grees of relativity and degeneracy), photons, and a mix- 
ture of non-relativistic classical Boltzmann gases of neu- 
trons, protons, a particles, and 14 kinds of heavy nu- 
clei. Coulomb lattice corrections of the pressure, energy 
density, entropy, and adiabatic index are taken into ac- 
count. Above a temperature of Tnse = 0.5 MeV the 
nuclear constituents are assumed to obey nuclear statis- 
tical equilibrium (NSE) , whereas below this temperature 
only nuclear burning of silicon and oxygen (and similarly 
neon and magnesium), and carbon (implemented accord- 
ing to the "flashing" treatment by Rampp & Janka 2002, 
Appendix B.2) can change the composition. The value 
of PeoS is chosen such that a smooth transition of pres- 
sure, internal energy density, and chemical potentials as 
functions of density is guaranteed. During the collapse 
phase peoS is set to 6 x lO^g cm~^ in case of the LS-EoS 
and 1.5 x lO^gcm"^ for the HW-EoS, while after core 
bounce, when a particles can reach a significant mass 
fraction and an error in the LS-EoS (Fryer, private com- 
munication) is potentially relevant (see, however, Buras 
et al. 2006a), we use pEoS = lO^^gcm"^. 

2.2. Investigated models 

Except for a simulation with an 11.2 progenitor 
from Woosley et al. (2002), for which we show results 
from a continuation of the model run published by Buras 
et al. (2006b), all core-collapse and supernova calcula- 
tions in this work are based on the 15Mq progenitor 
sl5s7b2 from Woosley & Weaver (1995) and are listed 
in Table 1. Details of the 11.2 M© model and a com- 
parison of ID and 2D results can be found in the Buras 
et al. (2006b) paper. Two of the 15 Mq calculations are 
conducted in spherical symmetry, one with the LS-EoS 
(indicated by "LS-ID" in the model name) and another 
one with the HW-EoS (labeled by -HW-ID"). These are 
compared with four 2D (axially symmetric) simulations, 
two of which (M15LS-2D and M15HW-2D) are computed 
without rotation, using the EoSs also employed in the ID 
models. In the other two simulations (i.e., our reference 
long-time run, M15LS-rot, and Model M15LS-rot9) we 



impose rotation on the progenitor core with initially a 
constant angular frequency of 0.5rads^^ (corresponding 
to a rotation period of about 12 s) inside the Fe core and 
an r~^/^ decline outside (see Fig. 1 in Miiller et al. 2004 
and Sect. 3.4 in Buras et al. 2006b). This rotation rate 
is roughly a factor ten faster than current predictions for 
pre-collapse stellar cores by Heger et al. (2005), but too 
slow for the free energy reservoir in the secularly evolv- 
ing proto-neutron star to be sufficient to power magneto- 
hydrodynamic explosions (see Burrows et al. 2007a and 
Thompson et al. 2005). 

In all of our nonrotating models we have implemented 
for the monopole of the gravitational potential the effec- 
tive relativistic potential according to Case A of Marek et 
al. (2006), which in that paper was identified to yield the 
best results in comparison to spherically symmetric, fully 
relativistic simulations. For models with core rotation, 
however, the best choice of the effective gravitational po- 
tential is less obvious. While for slow rotation Case A 
still produces the most satisfactory results, this monopole 
term is not able to account for the relativistic gravity of 
rotationally deformed cores. In such a situation, Marek 
et al. (2006) found Case R to be the better representa- 
tion, because, e.g., it leads to central densities closer to 
those of the relativistic calculation. Several 100 ms after 
core bounce even our slowly rotating pre-collapse model 
develops a sizable oblate deformation of the nascent neu- 
tron star. Since late-time effects are the main interest of 
our present investigations, we therefore decided to em- 
ploy the effective relativistic TOV potential of Case R of 
Marek et al. (2006) in our reference long-time 2D sim- 
ulation (Model M15LS-rot) and to compare this model 
with another simulation that includes rotation and uses 
Case A for the monopole term of the gravitational po- 
tential (Model M15LS-rot9). 

Our models are computed with 400-1000 non- 
equidistant radial zones, which are chosen such that the 
resolution Ar/r is typically better than three percent in 
the interior of the neutron star and between less than one 
percent and 1.5 percent around the neutrinospheres and 
outside of the neutron star. During the simulations the 
number of radial grid zones Nr is increased (Table 1 gives 
Nr at the beginning and at the end of the different runs) 
and the region of higher resolution is moved to smaller 
radii in order to follow the contraction of the neutron star 
and to ensure good resolution at the neutron star surface 
where the density gradient steepens with time. We make 
sure that at least 20 radial zones are used per density 
decade, corresponding to at most 12% density variation 
between neighboring zones. The representation of hydro- 
static equilibrium as well as the accurate computation of 
the neutrino transport in the decoupling layer require at 
least this numerical resolution; in case of the transport 
the luminosity variation from zone to zone is a limiting 
factor. Convergence tests for our code showed that insuf- 
ficient radial resolution in particular in this region may 
not only lead to bloated surface layers of the neutron 
star and to an overestimation of the neutrino luminosi- 
ties, but as a consequence also to artificial explosions 
even of spherically symmetric models. In order to make 
sure that our exploding 2D simulation. Model M15LS- 
rot, is well resolved, we recomputed a significant part 
of the crucial evolution phase of this model with about 
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half the radial zone width in the neutron star surface 
layer. This high-resolution run, Model M15LS-rot-hr, is 
not listed in Table 1, but some of the plots will also show 
results from this calculation for comparison. 

Our 2D simulations are performed with 128 or 192 
equidistant lateral zones from pole to pole of the spher- 
ical coordinate grid (Table 1). Comparing results for 
both resolutions we found excellent quantitative agree- 
ment of the neutron star and neutrino properties (which 
implies that the lateral resolution is sufficient for de- 
scribing neutron star convection) and very good qual- 
itative agreement of the shock evolution (because of 
the chaotic character of hydrodynamic instabilities in 
the hot-bubble region, a detailed quantitative agreement 
cannot be expected). At the grid center, however, a small 
region (1.7 km in radius or six radial zones) is computed 
in spherical symmetry to avoid the most severe limita- 
tion of the hydrodynamics timestep due to the Courant- 
Friedrich-Lewy condition. 

In the transport module of our code we usually employ 
17 geometrically distributed energy bins for neutrinos be- 
tween and 380 MeV. One of our models, M15LS-rot9, 
is computed only with nine energy bins, which leads to 
a reduced steepness of the high-energy tail of the neu- 
trino spectrum roughly three orders of magnitude and 
more below the spectral peak. Spherically symmetric 
test calciilations, however, show that this lower accuracy 
of the neutrino transport has hardly any influence on the 
hydrodynamic evolution of the supernova core and also 
causes only small differences (of order 10%) in the lu- 
minosities and mean energies of the radiated neutrinos. 
Using nine instead of 17 energy bins reduces the CPU 
time requirements of a simulation by nearly a factor of 
seven. 

The 2D runs without rotation are started at core 
bounce after mapping the corresponding ID models to 

the 2D grid and imposing random zone-to-zonc varia- 
tions of the radial velocity component (or density) with 
an amplitude of one percent. The rotating models are 
computed in 2D from the onset of the collapse with ran- 
dom perturbations and the initial rotation law imposed 
on the spherically symmetric progenitor star (for details, 
see Buras et al. 2006b). 

Running 2D models with a sophisticated energy- 
dependent treatment of the neutrino transport for many 
hundred milliseconds after bounce poses a considerable 
numerical challenge, because the rapid variations of the 
accretion flow in the presence of violent hydrodynamic 
instabilities in the postshock region do not allow the (im- 
plicit) neutrino transport timestep to become larger than 
about 2 X 10~^ s and therefore about 500.000 timesteps 
and 3 x 10^^ floating point operations were needed for 
our most advanced 2D model MS15LS-rot, which was 
evolved to more than 600 ms after bounce. For this rea- 
son, CPU-time limitations did not permit us to follow 
all of our 2D models to such very late post-bounce times. 
Although the set of calculations of Table 1 is not finished 
in this sense (and currently we have no way to do that 
in an acceptable time frame), it nevertheless allows for 
some interesting conclusions concerning the possibility of 
getting neutrino-driven explosions late after bounce and 
the dependence of this result on the properties and the 
neutrino emission of the forming neutron star. 



3. RESULTS 

In this section we discuss our simulations, starting with 

the ID results for reference, then continuing with the 
exploding long-time 2D simulation, and finally turning 
to the other 2D models, which allow us by comparisons 
to better understand the factors that are decisive for the 
explosion in the successful run. 

3.1. Spherically symmetric ISMq sim,ulations 

In Fig. 1 the results of our two ID simulations arc 
shown in overview. The collapse of the stellar core un- 
til bounce is somewhat faster with the softer (i.e. here: 
its adiabatic index in the subnuclear regime is lower) 
LS-EoS, and the maximum density reached at bounce 
is 3.6 X lO^'^gcm"^ compared to 3.0 x lO^'^gcm"^ in 
case of the HW-EoS. The central density is measured 
at a distance of 1.65 km (i.e., in the sixth radial grid 
zone) where numerical fluctuations caused by the inner 
grid boundary are absent). During the post-bounce evo- 
lution the higher stiffness of the nuclear phase in the 
HW-EoS leads to a considerably lower central density 
and a clearly less rapid contraction of the nascent neu- 
tron star, visible from the larger radius of the electron 
neutrinosphere in the right upper panel of Fig. 1. Since 
the postshock layer is nearly in hydrostatic equilibrium 
during the quasi-stationary accretion phase of the stalled 
shock, the shock position Rs is a sensitive function of the 
neutron star radius, Rns, roughly as 

/ 04 rr.4 

which follows from Eqs. (39), (56), and (63) of Janka 
(1991). Here M^s is the neutron star mass, T^s the tem- 
perature at the neutrinosphere, and M (< 0) the rate of 
mass accretion by the shock. In case of the HW-EoS, the 
stagnation radius of the shock is therefore significantly 
larger during the phase of shock retraction {t > 80- 
100ms after bounce). 

In contrast, the luminosities L,^ and mean energies 
(ejy) of the neutrinos radiated during shock accretion 
are appreciably higher in case of the LS-EoS, because 
the larger neutrinospheric temperature of the more com- 
pact neutron star overcompensates for the smaller radius 
(roughly, the emission behaves like blackbody radia- 
tion and thus oc R^^T^^ and (e^) oc T^s)- Interest- 
ingly, the prompt burst during shock breakout reveals 
the opposite dependence on the nuclear EoS: it is more 
luminous in case of the HW-EoS because of a stronger 
deleptonization in a wider spatial region, which is facili- 
tated by a less steep increase of the optical depth in the 
deleptonization region and thus an easier escape of the 
electron neutrinos. 

We also point out that in the simulations with both 
EoSs the average energy of the radiated z/g gets very 
close to that of the emitted muon and tau neutrinos or 
becomes even slightly higher after about 200 ms of post- 
bounce accretion (Fig. 1, bottom right). This effect is 
visible in the mean spectral energies, which are defined as 
ratio of the energy density to the number density of neu- 
trinos, (e^) = de Ji,(e)/ de e~ Vj.(e), with J^{e) 
being the zeroth energy moment of the specific inten- 
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sity"^. In contrast, the rms energies of the energy spec- 
trum, (e.)_ ^ [/~ de eV.(e)/ /~ de J.(e)] still fol- 
low the standard order sequence, < (epe)rms < 

(ej/ ) „, although the difference between the last two is 

considerably smaller than in older simulations, in which 
the transport treatment of heavy-lepton neutrinos Vx 
did not take into account the energy exchange through 
neutrino-nucleon scatterings and the production of u^Px 
pairs by nuclcon-nuclcon brcmsstrahlung and by the an- 
nihilation of z/ppp pairs (for more details, see Buras et 
al. 2003b, Raffelt 2001, Keil et al. 2003). We will come 
back to a closer discussion of these interesting spectral 
properties in Sect. 3.6. 

3.2. Two-dimensional 15Mq model with explosion 

In contrast to the spherically symmetric simulations, 
the two-dimensional Model M15LS-rot turns out to ap- 
proach an explosive runaway situation after more than 
500 ms of post-bounce accretion (Figs. 2 and 3). Some 
snapshots of the entropy distribution in the central region 
(with radii between ^400 km and ^800 km) for charac- 
teristic stages of the evolution are displayed in Fig. 4. 

3.2.1. Post-bounce evolution 

The results of a 90-degree equivalent of Model M15LS- 
rot (computed with a pole-to-equator wedge of the spher- 
ical coordinate grid and with reflecting conditions at both 
boundaries) until 300 ms after bounce were described by 
Buras et al. (2006b), see Sect. 3.4 there. Postshock con- 
vection is triggered by the small-scale perturbations that 
evolve from a 1% random initial (pre-coUapse) seed dur- 
ing infall. The convective activity exhibits the most 
rapid growth in high angular modes and visible inho- 
mogeneities begin to show up at 50-60 ms after bounce. 
In the 180-degree calculation discussed here, the first 
large-scale differences between both hemispheres appear 
at about 100 ms and only shortly later (at ~ 120 ms, see 
also the top left snapshot in Fig. 4) low-amplitude bipo- 
lar shock oscillations set in, manifesting the development 
of a global, low-mode asymmetry. This can be seen in the 
top panel of Fig. 5, where large amplitudes of the dipolar 
{I = 1) and quadrupolar {I = 2) terms of a spherical har- 
monics decomposition of the shock radius as a function of 
the polar angle 9 occur at i > 100 ms post bounce. The 
plot suggests that the dipole mode gains strength slightly 
faster and reaches large amplitudes some ten milliseconds 
earlier than the quadrupole mode. 

The appearance of strong low-mode, quasi-periodic 
shock deformation is associated with sloshing motions 
of the whole accretion layer between shock and neutron 
star. The process is very similar to the phenomenon iden- 
tified as SASI in the idealized accretion setup studied 
by Blondin et al. (2003) and seen also in more sophis- 
ticated post-bounce accretion simulations by Scheck et 
al. (2008). In those calculations convective activity was 
initially absent and set in only later, and therefore the 
oscillatory growth of the SASI in the linear regime could 

^ At sufBcicntly large radii, where neutrinos in the bulk of the 
spectrum propagate nearly radially (i.e., the flux factor in the lab- 
oratory frame is near unity, which is well fulfilled at the chosen 
radius of evaluation at 400 km), the local energy and number den- 
sities are essentially identical with the energy and number flux 
densities. 



be clearly identified. In contrast, in the model discussed 
here, sufficiently strong neutrino heating and sufficiently 
large seed perturbations cause convective activity in the 
gain layer to develop readily on a short timescale after 
bounce. The early presence of convection masks the char- 
acteristic growth behavior of the SASI. We suspect that 
the acoustic waves and vorticity created by convection 
even accelerate the evolution of the SASI and the condi- 
tions for a linear description never exactly apply. 

The growth of the SASI modes in the linear regime 
was analysed in detail by Blondin & Mezzacappa (2006). 
Their 2D axisymmetric hydrodynamic simulations with 
a simplified representation of the conditions in supernova 
cores revealed that the / = 1 mode is always unstable, 
grows faster than the I = 2 mode, and is the only mode 
that makes the transition to the nonlinear stage. Our 
simulations disagree with their results concerning the 
faster initial growth of the quadrupole mode (see the evo- 
lution during the first ^50 ms after bounce in the upper 
and middle panels of Fig. 5), which however was also seen 
in some model runs performed by Scheck et al. (2008). 
In the highly nonlinear phase we observe that the dipole 
and quadrupole modes arc present with similar strengths. 
During episodes in which the I = 2 component dominates 
the low-mode power (e.g. at t '--^150, 290, 380, 480 ms, 
see middle panel of Fig. 5), the spectrum is indeed en- 
hanced at frequencies well above 100 Hz (Fig. 5, lower 
panel), consistent with the expectation that the / = 2 
mode should be distinguished from the I = 1 mode by a 
higher frequency (Blondin & Mezzacappa 2006, Fig. 4). 
However, we cannot exclude that some of the activity 
we find associated with quadrupolar deformations is not 
a true I — 2 but a "ghost" I = 1 mode with the fre- 
quency of the latter, because the I = 2 spectra show a 
very broad distribution with significant power also below 
100 Hz. Perturbations associated with prompt postshock 
convection in the first tens of milliseconds after bounce 
(sec Marek et al. 2008) and violent neutrino-driven con- 
vection at t > 80 ms, which are present in all of our sim- 
ulations — different from the idealized accretion models 
considered by Blondin and Mezzacappa 2006, where the 
setup was chosen to be convectively stable — , but possi- 
bly also numerical errors must be expected to give rise to 
the excitation of Z = 2 motions and thus might prevent 
the presence of a purer and more clearly dominant I = 1 
mode. 

In Fig. 3 the the radii of the shock near the north pole 
and south pole are plotted as functions of time. One can 
well sec the alternating phases of shock expansion in the 
northern and southern hemispheres. The snapshots in 
Fig. 4 at 454 ms and 524 ms after bounce capture two 
such moments in which the shock is inflated in opposite 
directions. There is an inverse relation between the typi- 
cal period of these SASI oscillations and the mean shock 
radius. The period is about 10-12 ms when the shock is 
near 150 km (around 300-400 ms after bounce) and more 
like 15-17 ms when the average shock radius reaches 180- 
200 km (around 500 ms), corresponding to power maxima 
at typical frequencies beween 50 Hz and 100 Hz. Towards 
the end of the simulation, as the shock expands farther, 
the oscillation period becomes even longer. This can be 
seen in the time-dependent frequency spectrum of the 
low-mode SASI power (obtained from a fourier analysis 
of the quantity given in Eq. (2)) in the bottom panel of 
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Fig. 5. The power spectrum has a very broad peak in the 
mentioned frequency interval and decays steeply towards 
higher frequencies. 

The middle panel of Fig. 5 displays the corresponding 
time evolution of the total power in the I = 1 and I = 2 
modes, computed as the volume integral of the squared 
dipole and quadrupolc amplitudes of the fractional pros- 
sure variations between average neutrinosphere radius 
Rv and average shock radius Rg (cf. Blondin & Mez- 
zacappa 2006), 

Power(Z;i) = 27r / dr r'^ of o{r,t) , (2) 

where aifi{r,t) are the amplitudes of the spherical har- 
monics expansion of the normalized pressure fluctuations 
according to 

= f (3) 

\ y ^ T )/0 

with ^^(cos^) being the Legendre polynomials (the in- 
dex for the azimuthal modes is m = because of the 
axial symmetry of our 2D models). The combined power 
of the two lowest modes turns out to reflect rather sensi- 
tively the dynamical activity in the accretion layer. One 
should note that during phases of relative quiescence of 
the dipole mode the quadrupole mode is dominant and 
vice versa. A high level of activity is reached shortly 
after convection has become strong and the SASI defor- 
mation of the shock has set in {t > 100 ms after bounce). 
The following slight reduction of the power is a conse- 
quence of the shock retraction between 100 and 150 ms 
post bounce. When the jump in the entropy, density, 
and mass accretion rate associated with the composition 
interface between the Si-layer and the oxygen enriched 
Si-shell of the progenitor reaches the shock at 170 ms 
after bounce, transient shock inflation is triggered (cf. 
Figs. 1, 2, 3, and 6). As a consequence, the SASI power 
increases again before it decays once more during an- 
other period of shock contraction. At f > 400 ms a phase 
of basically continuous, slow expansion of the average 
shock radius begins and the low-mode power grows. Af- 
ter 500 ms until the end of our simulation at ^700 ms, 
the low SASI modes gain another factor of 100 in power 
(middle panel of Fig. 5), which is accompanied by an ac- 
celerating mean shock inflation. Finally, after 550 ms, the 
expansion of the average shock radius speeds up strongly, 
and Model M15LS-rot enters a runaway situation that 
initiates an explosion (Figs. 2, 3, 4, and 6). The shock 
acceleration is supported by the onset of nucleon recom- 
bination to a-particles in the hot-bubble medium or/and 
by only partial dissociation of a-particles in swept-up 
matter as the shock reaches radii beyond 200 km and 
temperatures of less than ~1 MeV are present behind 
the shock (cf. the light-grey region behind the shock in 
Fig. 2). A nuclear binding energy of about 7 MeV per 
nucleon is either released by the formation of a-particles 
or saved when helium is not dissociated. This increases 
the thermal energy of the stellar gas and thus adds to 
the energy deposited by neutrino heating in bringing the 
gas in the gain layer to a gravitationally unbound state 
(for a further discussion of this point, see Sect. 3.4). 

During the outward acceleration of the average shock 
radius at t > 550 ms, the two-dimensional shock con- 
tour in the axisymmetric Model M15LS-rot develops a 



pronounced butterfly-like shape (middle right and bot- 
tom panels of Fig. 4) with alternatingly stronger expan- 
sion in either the northern or the southern hemisphere 
(Fig. 3) and a significantly more infiated southern lobe 
at the end of the simulated evolution. The nearly unipo- 
lar onset of the explosion indicates the dominance of the 
I = 1 SASI mode, consistent with the fact that the am- 
plitude of the I — 1 component is clearly dominant in 
the upper and middle panels of Fig. 5 at f > 600 ms. 
The swelling high-entropy plumes of neutrino-heated gas 
in both hemispheres push the shock outward and cre- 
ate a kink in the shock surface near the equator. The 
kink feature causes accreted gas to hit the shock at an 
oblique angle and therefore to be deflected at passing 
the shock. The deflected gas streams from the northern 
and southern hemispheres collide and form a long-lived 
and waving equatorial downdraft. This accretion funnel 
collects the major fraction of the equatorially infalling 
gas and channels it towards the neutron star with su- 
personic velocities. The existence and stability of this 
feature has important consequences for the long-lasting 
energy input by neutrino heating into the developing ex- 
plosion (see Sect. 3.4). Gas with low angular momentum, 
however, collapses closer to the rotation axis and creates 
very stable polar downflows in both hemispheres. These 
polar downdrafts are a common feature of all 2D simula- 
tions with rotation and their durability and universality 
can be understood by the Solberg condition, which says 
that a situation in which the specific angular momentum 
increases with distance from the rotation axis is stable 
(see the results and discussion in Scheck et al. 2006). 

Rotation, although "modest" in the sense discussed in 
Sect. 2.2, has an important influence in Model M15LS- 
rot, but it is not crucial for the evolution after bounce and 
the final outcome obtained in this simulation. This will 
be further discussed in Sect. 3.5 by comparing the diff'er- 
ent models of our set of runs. Until the end of the compu- 
tation (at ~615ms), the neutron star in Model M15LS- 
rot has contracted to an average radius of 26 km (com- 
pared to 21 km for the ID Model M15LS-1D) with a pole- 
to-equator radius ratio of 0.65 (eccentricity of 0.76). It 
rotates strongly differentially with an average period of 
3 ms. In the gain layer the rotation period is typically 
longer than 10 ms. 

3.2.2. Runaway conditions and explosion criteria 

At t > 350 ms Model M15LS-rot evolves gradually and 
steadily towards a situation that is more and more favor- 
able for an explosion. This is clearly visible in a number 
of quantities. Not only the continuous expansion of the 
shock front (Figs. 2, 3, and 6) but also the mass, en- 
ergy, and neutrino heating in the gain layer reveal the 
trend that the explosion conditions become better at 
later times. 

Two timescales are crucial for the behavior of the gain 
layer: the timescale of accreted matter to be advected 
from the shock to the gain radius. 



and the neutrino heating timescale of the gas in this re- 
gion, 

_ |^^bind|[-Rg, Rs] /rx 

Theat = (5) 

*^heat 
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Here Rg denotes the average gain radius, (vr) is the av- 
erage postshock velocity, ii'bind the total energy of the 
matter in the gain layer as sum of internal, kinetic, and 
gravitational energies, and Qheat is the integrated net 
heating rate by neutrinos. Therefore Theat measures the 
time it takes neutrinos to deposit an energy equal to the 
binding energy -Ebind [-Rgain , -Rshock] of the matter in the 
gain layer, whereas Tadv can be considered as a measure 
of the duration the gas is exposed to neutrino energy 
deposition. In multi-dimensional simulations when non- 
radial hydrodynamic instabilities play an important role, 
Eq. (4) for evaluating the advcction timescale should be 
replaced by Eq. (8) of Buras ct al. (2006b). A ratio 
of Tadv/Thcat > 1 Can bc considcrcd as indicative of fa- 
vorable conditions for an explosion (Janka et al. 2001, 
Thompson et al. 2005, Buras et al. 2006b). Of course, the 
definitions of both timoscalcs contain a variety of ambi- 
guities and uncertainties, e.g. whether the total energy or 
just the internal energy is the more suitable energy scale 
in the nominator of Eq. (5) or where and how exactly the 
mean advection velocity in the denominator of Eq. (4) 
should be measured. Therefore the timescale ratio is 
only a rough diagnostic number that can be indicative 
of the presence or absence of a generally favorable trend, 
but certainly Tadv/Theat > 1 cannot be taken as a rigor- 
ous quantitative criterion for the onset of the successful 
blast. Murphy & Burrows (2008) also stress this fact and 
instead recommend to monitor the residence time distri- 
bution function of tracer particles, which is better able to 
capture the complex flow dynamics in multi-dimensional 
simulations. While this is unquestionably correct, the 
numerical effort of tracking tracer particle paths is con- 
siderably bigger than performing an analysis of the simu- 
lations on the basis of Eqs. (4) (or its multi-dimensional 
analog as given by Eq. (8) of Buras et al. 2006b) and 
Eq.(5). 

Figure 7 shows that along with the growing strength of 
the low SASI modes (cf. Fig. 5) the advection timescale 
in Model M15LS-rot rises continuously after 350 ms post 
bounce, leading to a steady growth of the timescale ra- 
tio Tadv/Theat from about 0.6 to values larger that unity. 
The critical value of unity is exceeded at '^530ms after 
bounce and the ratio remains larger than unity until the 
end of the simulation, i.e. for a period of time much 
longer than the heating timescale Theat ~ 20 ms. This is 
a necessary condition for an explosion (see Buras et al. 
2006b). Note that even after the shock starts a rapid ac- 
celeration at 550 ms, Tadv can be evaluated because there 
is net accretion still going on. At the end of our simu- 
lation for Model M15LS-rot, about 50% of the matter 
swept up by the shock is still accreted onto the neutron 
star (the fraction is larger at earlier times), while the 
rest stays in the gain layer and causes a growth of the 
mass there at a rate of more than 0.1 AIq s~^ (see the top 
panel in Fig. 8). The dominance of mass inflow towards 
the neutron star compared to mass outflow through the 
gain radius until t > 650 ms after bounce cannot only be 
concluded from the fact that the advection timescale is 
still defined at such a late time. It is also suggested by 
the mass shell trajectories in Fig. 2 at this time, which 
indicate that a significant fraction of the matter in the 
gain layer is still crossing the gain radius with negative 
velocities. 

The neutrino-heating timescale of the gain layer in 



Model M15LS-rot levels off near 20 ms after t ^ 350 ms 
(Fig. 7, middle panel). This corresponds to a net rate of 
neutrino energy deposition (the "net rate" is defined as 
the rate of energy input by neutrinos minus the rate of en- 
ergy loss by neutrino reemission) of about 4 x 10^^ erg s~^ 
(Fig. 8, middle panel) with a mild increase with time. 
This increase is a conscqiience of the slowly growing mass 
and optical depth of the gain layer and of the higher 
mean energies of the neutrinos radiated at later times 
(see Fig. 1 and 14). The corresponding neutrino heating 
efficiency, defined as ratio of the total net energy depo- 
sition to the sum of the radiated ly^ and 9^ luminosities, 
rises from about 6% at 350 ms to peak values of more 
than 8% at times later than 500 ms after bounce. 

With the accumulation of mass in the gain layer, also 
the internal energy in this region grows (Fig. 9, upper 
right panel). The same is true for the rotation energy of 
the accreted matter, which rises becaiisc the gas that is 
accreted at later times comes from larger distances and 
has a higher specific angular momentum. The increase of 
both energies, however, is rather moderate (some 10%) 
compared to the growth of the kinetic energy of fluid mo- 
tions in the gain layer in non-azimuthal directions. Both 
the energy associated with gas flows in latitudinal direc- 
tion and with gas motions in radial plus lateral directions 
gain a factor of roughly four between 350 ms and 600 ms 
(Fig. 9, upper left plot). This, of course, reflects the in- 
creasingly more violent SASI and convective activity in 
the shocked flow. These energies are around 10^'' erg, 
corresponding to ~0.7MeV per nucleon, a value that is 
dwarfed by the neutrino heating rate (300MeVs~^ per 
nucleon), applied for a typical timescale Tadv, and by 
the recombination energy of nucleons to a-particles. Ac- 
cordingly, the total energy per nucleon in the gain layer 
changes much more strongly from a mean value of about 
-7MeV at 350 ms to -3.5 MeV around 600 ms and to 
only -2 MeV at 700 ms. 

Towards the end of the computed evolution of 
Model M15LS-rot a sizable and steeply rising fraction 
of the mass in the gain layer has obtained a specific en- 
ergy that is either positive or only marginally negative 
(Fig. 9, lower two panels). The energy of all mass in the 
gain layer with positive total (internal plus kinetic plus 
gravitational) specific energy and positive radial velocity 
at that time is roughly 2.5 x 10^^ erg and also steeply 
increasing (Fig. 9), correlated with a very steep rise of 
the kinetic energies near the termination point of the 
simiilation. 

A variety of energy parameters therefore signals that 
an explosion begins to develop in Model M15LS-rot. This 
can also be concluded from the behavior of the advection 
timescale and the timescale ratio at the end of our sim- 
ulation. Both show a steep increase because more and 
more matter in the gain region attains positive velocities 
and starts expanding (Fig. 7). Since a growing fraction 
of the matter that is accreted by the shock stays in the 
gain layer, the mass in this layer grows and the gas infall 
rate onto the neutron star is reduced. Correspondingly, 
less neutrinos are radiated from the accretion layer near 
the neutron star surface. This leads to a visible reduction 
of the luminosities of Ve and De and of the mean energies 
of all kinds of neutrinos (Fig. 14). 

3.2.3. Role of the SASI 
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In the last section we argued that the combination 
of shock expansion, critical timcscalc ratio, and grow- 
ing mass and energy in the gain layer indicates that 
Model M15LS-rot approaches a runaway situation. The 
described behavior of these quantities is very similar to 
the findings of Scheck et al. (2008) for their exploding 
model WOOF (cf. Figs. 19 and 20 in their paper). In the 
latter model the growth of buoyancy instabilities was ini- 
tially damped due to the choice of small seed perturba- 
tions and a low neutrino-heating rate in the gain layer. In 
contrast, the larger perturbations that have developed in 
the stellar core of Model M15LS-rot during the collapse 
allow convection to begin in this simulation as soon as 
neutrino heating achieves to create a convcctively unsta- 
ble situation in the gain layer. Therefore the growth of 
the SASI and of convection cannot be studied separately 
in the present simulation. 

Nevertheless, the large-amplitude low-mode shock os- 
cillations seen in Model M15LS-rot only ^^lOOms after 
bounce (see Figs. 3 and 5) clearly indicate the presence 
of the SASI, whose growth appears to be accelerated by 
the sonic and vorticity perturbations associated with the 
convective activity in the postshock region. Systematic 
studies by Scheck et al. (2008) showed that the SASI can 
play an important supportive role for neutrino-driven ex- 
plosions. Scheck et al. (2008) demonstrated that SASI 
shock oscillations can develop even in conditions which 
per sc are disfavorable for the growth of convection. Thus 
the SASI can assist the delayed revival of the stalled 
shock in at least two ways: (i) Large-amplitude expan- 
sion and contraction phases of the shock lead to super- 
sonic lateral velocities in the postshock flow and the for- 
mation of sheets with very steep unstable entropy gra- 
dients. This fosters the development of secondary con- 
vection or strengthens the ongoing convective activity, 
(ii) SASI sloshing motions as well as strong postshock 
convection do not only produce large non-radial veloc- 
ity components in the postshock flow but also push the 
accretion shock to larger radii. This again reduces the 
average infall velocity in the postshock layer, because the 
velocity in the free-fall region ahead of the shock drops 

— 1/2 

like Rs .A larger shock radius thus leads to a signifi- 
cantly longer advection timescale, roughly Tadv oc Rs^^- 

These effects help to keep accreted matter in the gain 
layer for a longer period of time (a consequence of which 
is the increasing mass in this layer). This supports the 
energy deposition by neutrinos and thus facilitates the 
explosion. So, although the kinetic energy associated 
with the SASI remains negligible for the explosion en- 
ergetics, once the SASI and convection have pushed the 
shock sufficiently far out, prolonged neutrino heating and 
ultimately also the additional energy release by nucleon 
recombination (or, similarly, energy savings by only par- 
tial dissociation of nuclei in the shock-accreted gas) set 
up the conditions for runaway. 

The contributing effects in the postshock accretion 
flow, the SASI, convection, and neutrino heating, how- 
ever, are strongly linked and interdependent in a complex 
way. In such a highly nonlinear situation it is therefore 
extremely difficult to exactly determine the influence of 
each aspect individually. It is, for example, unclear which 
conditions are primarily responsible for the growth of the 
SASI power at late post-bounce times {t > 300 ms), and 



which role neutrino heating may play in this context. By 
comparing in Sect. 3.5 Model M15LS-rot with the other 
simulations of our sample — although it was not yet 
possible to carry these other runs to equally late times 
as our reference model — we will strive for a better un- 
derstanding of the influence of at least some potentially 
important aspects like rotation and the contraction and 
growing compactness of the forming neutron star. 

3.3. Two-dimensional 11.2Mq m,odel with explosion 

As discussed in the previous sections, the strong am- 
plification of the SASI and convective activity lead to 
runaway conditions and the initiation of an explosion 
in Model M15LS-rot for a ISM© star. Buras et al. 
(2G06b) observed a similar situation in their 180-degree 
simulation of the collapse and post-bounce evolution of 
a (non-rotating) 11.2 Mq progenitor from Woosley et al. 
(2002). Different from our successful 15 Mq run, the 
neutrino-powered explosion was found to set in much ear- 
lier, namely already at about 200 ms after bounce. 

The importance of the presence of low-mode SASI os- 
cillations for the development of the blast in this case 
was concluded from the comparison with a simulation 
that was computed with the same input physics but with 
the lateral grid constrained to a 90-degree wedge around 
the equatorial plane of the polar coordinate system and 
periodic boundary conditions. This simulation did not 
show an explosion until the same post-bounce time''. In 
the latter setup the choice of the grid did not allow for 
the occurrence of / = 1 and I = 2 SASI modes. In- 
terestingly, Marek (2007) obtained an explosion for the 
11.2 Mq progenitor when he reran the model with a 90- 
degree grid extending from the pole to the equator with 
reflecting boundaries, thus assuming equatorial symme- 
try and excluding the I = 1 mode but giving room for the 
mirror symmetric / = 2 mode. This result was recently 
confirmed by Murphy & Burrows (2008) and shows that 
obviously the support of strong quadrupolar SASI mo- 
tions makes already a significant difference compared to 
a case where both lowest-order spherical harmonics are 
suppressed. 

Here we present some results for a continuation of the 
11.2 M© explosion simulation of Buras et al. (2006b). 
Due to the larger timesteps that were possible in this 
case, we could evolve the model further into the explo- 
sion. The results are therefore of interest in the present 
context because they provide insight into the way the 
explosion strengthens at evolution stages that we were 
unable to reach in the 15 Mq run. 

Figure 10 displays the mass-shell evolution as deduced 
from the (mass- weighted) angle-averaged 2D data of the 
simulation. Grey or yellow shading of different regions 
gives rough information about the nuclear composition 
in the corresponding layers. One can see that the onset 
of the shock expansion coincides with only partial disso- 
ciaton of helium nuclei to free nucleons in the postshock 
region. Whether this is causal for the shock expansion or 
just a consequence of it, is not finally clear. In any case, 

^ Of course, this does not permit the conclusion that an ex- 
plosion could not occur later. Such a possibility, however, seems 
disfavored because the sensible parameters for a runaway situation 
(advection timcscalc, neutrino heating rate, mass in the gain layer, 
etc.) exhibit a much more pessimistic trend than in the successful 
11.2 Mq run with the 180-degree grid. 
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the larger abundance of a-particles reduces the consump- 
tion of thermal energy for breaking up strongly bound 
nuclei. Figure 10 also shows that the shock reaches the 
infall layers containing a mass fraction of more than 10% 
oxygen at f « 270 ms after bounce only after it has al- 
ready propagated to more than 700 km. 

At that time the shock is already half way through the 
supersonically collapsing inner shells of the progenitor 
star and on average rushes outward with a radial velocity 
of 10,000 kms~^ (Fig. 11, left panel). The shock surface 
has a highly deformed prolate shape with an z : x-axis ra- 
tio of more than 2:1 (Figs. 11 and 12). Large lobes filled 
with neutrino-heated, high-entropy gas expand into the 
northern and southern hemispheres and arc continuously 
fed with fresh material that is channelled in a long-lived, 
very stable equatorial accretion funnel from behind the 
shock to the close vicinity of the gain radius to absorb en- 
ergy there from neutrinos. The situation is very similar 
to what we observed in the case of our 15Mq explosion 
model (see Sect. 3.2.1). Until the end of the simulated 
evolution, the outward moving, hot gas in the two big 
lobes (the northern one being slightly larger) has accu- 
mulated a positive total (internal plus kinetic plus gravi- 
tational) energy of nearly 2.5 x 10**^ erg with a tendency 
of steep rise (Fig. 11, right panel). 

The pronounced asphericity has the consequence that 
the shock reaches the silicon-oxygen interface near the 
equatorial plane later than in the polar directions (see the 
panels for t = 250 ms and 275 ms after bounce in Fig. 12). 
Therefore a wedge-like region around the equator re- 
mains for some time, where silicon and sulfur are still 
present with higher abundances between the shock and 
the oxygen layer, while the matter swept up by the shock 
consists mostly of iron-group nuclei and a-particlcs. The 
mass-shell plot of Fig. 10, which is constructed from the 
laterally averaged 2D data at each radius, is mislead- 
ing by the fact that this preshock material appears to 
be located behind the angle-averaged shock radius (at 
post-bounce times 270 ms ^t< 300 ms). We note that 
the penetration into the oxygen-rich infalling shells, be- 
ginning at t ~ 250 ms p.b., does not have any obvious 
supportive or strengthening effect on the outgoing shock. 

In Fig. 13 we provide information about the condi- 
tions and neutrino-energy deposition in the gain layer 
of the 11.2 M0 model. As in the 15 Mq case the mass 
in the gain layer increases when the shock begins its 
outward expansion. At the same time the infall (ad- 
vection) timcscale of matter between the shock and the 
gain radius increases, but continues to be well defined. 
Again, as in the 15 Mq explosion model, this suggests 
the presence of ongoing accretion of gas through the 
gain layer to the neutron star (which can also be con- 
cluded from the continued contraction of mass shells in 
this region in Fig. 10). Shortly after the (net) neutrino- 
heating rate has reached a pronounced peak of about 
7.5 X 10^^ ergs^^ at i w 70ms, it makes a rapid drop to 
around 3 x 10^^ ergs~^. This decline is a consequence of 
the decay of the neutrino luminosities at the time when 
the mass infall rate onto the shock and the neutron star 
decreases. The decrease occurs when the steep nega- 
tive density gradient (and positive entropy step) near 
the composition interface between the silicon layer and 
the oxygen-enriched Si- layer of the progenitor star (near 
1.3 Mq) arrives at the shock (at t w 100 ms after bounce). 



Nevertheless, the heating timescale shrinks essentially 
monotonically, which points to an evolution of the mat- 
ter in the gain layer towards an unbound state, i.e., the 
absolute value of the total gas energy in the numerator 
of Eq. (5) goes to zero. 

3.4. Explosion energy 

In both our 11.2 and 15 Mq explosions, the en- 
ergy of the matter in the gain layer with positive radial 
velocities ("explosion energy") reaches '^2.5x 10^^ erg at 
the end of the computed evolutions and rises with a very 
steep gradient (Figs. 9 and 11). Therefore reliable esti- 
mates of the final explosion energy cannot be given at 
this time. For that to be possible, the simulations would 
have to be continued for many himdred milliseconds more 
(which is numerically a challenging task and currently 
impossible for us with the sophisticated and computa- 
tionally expensive neutrino transport and chosen resolu- 
tion). This is obvious from the neutrino-driven explosion 
models investigated by Scheck et al. (2008; Figs. 9 and 
10 there) and Scheck et al. (2006; Appendix C there) 
and has several reasons. (1) When rapid shock acceler- 
ation sets in, only a smaller fraction of the mass in the 
gain layer starts to attain positive total specific energy 
(see Fig. 9 for our 15 Mq explosion model and Fig. 24 
in Buras et al. 2006b for the 11.2 simulation), and it 
take il ten milliseconds until a major part of the 

mass between shock and gain radius follows the outward 
motion of the shock and the explosion energy reaches 
even only some 10*® erg (cf. Fig. 10 in Scheck et al. 2008, 
Fig. 11 and Fig. 24 in Buras et al. 2006b for the 11.2 Mq 
explosion, and Fig. 9 for the exploding 15 Mq model). 
(2) Also then the energy of the explosion grows only 
gradually and initially roughly linearly with a rate of 
some 10^° erg per 100 ms. During this phase ongoing 
accretion transports fresh gas towards the gain radius, 
where a part of the gas absorbs energy from neutrinos 
and begins to rise again, while the rest of the accreted 
gas is advected into the cooling layer below the gain ra- 
dius and is added into the forming neutron star. (3) 
After some time the downdrafts of accreted matter may 
be quenched and the inflow of gas towards regions near 
the gain radius may stop. Nevertheless, more neutrino- 
heated gas is ejected in the neutrino-driven wind that 
sheds off matter from the surface layer of the nascent 
neutron star. During this phase the increase of the ex- 
plosion energy levels off, but still a signiflcant fraction of 
the final energy may be added (though with a much lower 
rate) over timescales of several more seconds. The corre- 
sponding wind power is a very sensitive function of the 
neutrino emission properties (i.e., of the time-dependent 
luminosities and mean energies of the radiated neutrinos) 
and of the neutron star mass and radius. Some 10^" erg, 
in optimistic estimates even more than 10^^ erg, of en- 
ergy might be pumped into the supernova this way after 
the onset of the blast (see Burrows & Goshy 1993; Qian 
& Woosley 1996; Thompson, Burrows, & Meyer 2001). 
(4) The total energy injected at the explosion center has 
to be corrected for the negative total energy of the grav- 
itationally bound outer stellar shells, which are going to 
be swept out in the blast instead of falling back to the 
neutron star. 

Point (2) implies that simultaneous accretion and 
outflow of neutrino-heated g clS clS cl generically multi- 
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dimensional phenomenon is an essential feature of the 
supernova explosion mechanism. The rate at which 
accreted, cool gas is channelled through the neutrino- 
heating region, gains energy from neutrinos, and finally 
rises outward again (the process described in point (2) 
above) decides about the power that is accumulated in 
ejected material. Let the corresponding outflow rate of 
neutrino-heated gas be a fraction ( of the mass accretion 
rate Mace through the shock, i.e., Mout = C-^acc- The 
rate at which energy is pumped into the ejecta can then 
be estimated very roughly as 

~ C-^acc^iyTadv 

^""^^ s [o.5j \0.2Mq/s) 

^ 1^300 MeV/sy V30ms/ ' ^ ^ 

Here tub is the baryon mass and q^, the average net neu- 
trino energy deposition rate per unit of mass in the gain 
layer, which can be estimated by the ratio of the net 
neutrino heating rate in the gain layer, E'l^^gain (middle 
panel of Fig. 8), to the mass Mgain in the gain layer (up- 
per panel of Fig. 8): 

^.«|^. (7) 

The advection timescale Tadv through the gain region in 
Eq. (6) is a measure of the mean exposure time of the 
gas to neutrino heating^. The numerical values used in 
Eq. (6) are guided by the conditions a.t t > 600 ms in 
run M15LS-rot (see iPigs. 6, 7, 8, and Sect. 3.2.2), for 
which we therefore estimate Ei, ~ (1--2) x lO^^ergs"^. 
In the case of our explosion simulation of the 11.2 M© 
star, Eq. (6) with the numbers taken from the latest evo- 
lution stage shown in Fig. 13 yields a neutrino energy in- 
put rate to the ejecta of E^, ~ (3. ..6) x 10^° ergs~^, which 
is consistent with the growth of the positive energy in the 
gain layer at the end of the model run (compare Fig. 11, 
right panel). The energy deposition rate according to 
Eq. (6) is of the order of the net neutrino heating rate of 
the gain layer, £^i,,gain (see Figs. 8 and 13). In contrast 
to the latter, howcn^er, it accounts for the fact that not 
all of the matter that absorbs energy from neutrinos in 
the gain layer will finally be added to the outflow (see 
also Sect. 3.2.2). We note that the power of the outflow 
of neutrino-heated matter is seriously underestimated in 
simulations with Newtonian instead of relativistic grav- 
ity, because in the Newtonian case the luminosities and 
energies radiated during the accretion phase of the less 
compact and cooler neutron star are considerably lower 
and Cfi, is correspondingly smaller. 

With E^ from Eq. (6) and the duration Tacc of the 
accretion phase, the explosion energy, in a crude way, 
can be written as 

-^exp ^ -^i/'^acc ~t~ -^wind ~t~ -^burn -^bind * (S) 

^ The product Mace Tadv turns out to be a fairly accurate rep- 
resentation of the mass in the gain layer at all times after bounce. 
This fact confirms that the advection time according to the defini- 
tion in Eq. (8) of Buras et al. (2006b), which is used to evaluate 
Tadv ™ our models, is reasonable and meaningful. 



Here the first term on the rhs side accounts for the to- 
tal energy that is transferred by neutrinos to the ejecta 
during the phase of simultaneous accretion and out- 
flow. The second term measures the integrated power 
of the neutrino-driven wind that starts after accretion 
has ceased (the phase mentioned in point (3) above). 
Wind matter in contrast to accreted gas must be lifted 
all the way out of the deep gravitational potential at the 
neutron star surface. The third term takes into account 
additional energy release by nuclear burning in shock- 
heated matter when the shock wave reaches layers of still 
unburned nuclear fuel, and the fourth term corrects for 
the gravitational binding energy of the overlying stellar 
material (as addressed in point (4) above). For stars 
around 15-20 Mq the last two terms can be estimated to 
approximately compensate each other (see Janka et al. 
2001). In the neutrino- driven wind phase (when accre- 
tion of infalling gas onto the neutron star has ended) the 
gain radius retreats to a location closely outside of the 
neutrinospheres (this phase is not yet reached in both the 
11.2 Mq and 15 Mq explosions discussed in this paper, 
see Figs. 2 and 10). 

The release of nuclear recombination energy E^ec when 
free nucleons in the neutrino-driven wind ejecta assem- 
ble to a-particles and heavier nuclei, must be included 
as a positive contribution to the £^wind term so that 
^/wind = -E^.gain + E^cc during the wind phase. For mat- 
ter that is accreted through the shock and gets heated 
by neutrinos (an effect that is accounted for in the first 
term on the rhs of Eq. (8)), however, such nucleon recom- 
bination does not necessarily mean a gain in the energy 
budget. The gas swept up by the shock and falling inward 
to the neutron star in downflows is initially composed of 
nuclei (iron, silicon, or oxygen, depending on the location 
of the shock, see Fig. 2). Compressional heating within 
the shock and in the downflows (or ultimately energy de- 
position by neutrinos) produces such high temperatures 
that the nuclei get disintegrated, which consumes a fair 
amount of dissociation energy. This energy can subse- 
quently be recovered in the fraction of the matter that 
reexpands and gets blown out again. Different from the 
energy input by neutrinos, the recombination of nucle- 
ons, however, yields a net energy gain only when the 
nucleon recombination leads to a state with higher nu- 
clear binding energy per baryon. This means that the 
net gain of energy by nuclear processes in the accreted 
material can at most be equal to the difference in the 
binding energy of oxygen or silicon (initially present in 
the accretion flow) compared to iron-group elements (the 
most strongly bound nuclei), which is much smaller than 
the complete recombination energy of free nucleons to 
nuclei. 

For both components of the ejecta, however, for the 
neutrino-driven wind material as well as accreted and 
reexpanding matter, nuclear processes are unlikely to be 
the main source of the energy excess (and thus of the ex- 
plosion energy) that the matter carries outward in the su- 
pernova blast. On the one hand, gas residing in the gain 
layer before the explosion sets in as well as matter ablated 
from the neutron star surface in the later neutrino-wind 
phase is initially at rest in the strong gravitational field 
very close to the neutron star and has a largely negative 
total energy (©(— 10 MeV/nucleon) or more). The en- 
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ergy liberated by the recombination of free nucleons to 
nuclei in this matter, though large (^7 8.8MeV per nu- 
clen), can typically only assist the revival of the stalled 
shock but cannot account for all the energy needed to 
make this matter gravitationally unbound. Neutrino en- 
ergy deposition is also needed to overcome the gravi- 
tational binding and to produce the positive energy of 
the ejecta at infinity. On the other hand, gas that is 
accreted from large distances, channelled through the 
neutrino-heating region, and then reexpanding outward 
right away is initially gravitationally only weakly bound 
(its total specific energy is typically only of the order of 
— (1...2) MeV per nucleon). A change of the nuclear com- 
position in this matter by assembling heavier nuclei from 
lighter ones can release only a relatively small amount of 
energy (also of the order of 1-2 MeV per nucleon) and can 
bring the gas at most to a gravitationally marginally un- 
bound energy level. So again the ultimate energy excess 
that drives the supernova explosion has to come from 
another source. Neutrino energy deposition plays this 
dominant role in our simulations. 

In this context it is crucial to answer the question how 
long accretion can continue and how long efhcient energy 
input by neutrinos according to Eq. (6) can therefore 
apply. The time when accretion ends can be deduced in a 
crude way from the requirement that at this moment the 
outgoing shock (with velocity Vs and postshock velocity 
V2 in the observer frame) must accelerate the postshock 
matter to escape velocity. The corresponding radius is 
therefore determined by 

2GM o 
Resc ^ — — ^ 5400km Mi,5V-^ , (9) 

where Mi,^ denotes the enclosed mass normalized to 
1.5Mq, ?;s_9 the shock velocity in lO^cms^^, and the 
postshock velocity V2 was determined from the preshock 
velocity Vi and the shock velocity Vs by the first shock- 
jump condition as 

V2 = ii-r')vs + r'vi ^ ^vs, (10) 

with /3 = P2/P1 ~ 7 and vi ~ 0. For the duration of the 
accretion phase after the onset of outward shock expan- 
sion we therefore get: 

race « — ~ 0.5sMi.5V-9^ . (11) 

This means that accretion can easily continue for half a 
second or even longer, if the shock expands slowly in the 
beginning. This is confirmed by the supernova simula- 
tions of Scheck et al. (2006), in some of which accretion 
was observed to still go on significantly longer than one 
second after the explosion had started. The explosion en- 
ergy is therefore not determined by the conditions just at 
the moment when the outward shock acceleration begins 
or when the shock is in its early stage of expansion. For 
this reason it is impossible to conclude from our simula- 
tions that the 11.2 Mq and 15 Mq explosions will finally 
be weak or strong. In the multi-dimensional environment 
of the supernova core with accretion and shock expan- 
sion being simultaneously present, it is essential to take 
into account the long-lasting period of mass downflow to 



the nascent neutron star^. Even though the neutrino- 
driven explosion sets in rather late after bounce in our 
simulation, Eqs. (6) and (11) suggest that the nonspher- 
ical nature of the gas flow in the supernova core, due 
to which the neutrino-heating region around the neutron 
star is replenished with fresh gas over a long period of 
time, may well allow the explosion energy E^xp ~ E^t^cc 
to reach several 10^° erg for the 11.2 Mq progenitor and 
even a value around the canonical number of 10^^ erg in 
the case of our 15 Mq explosion model. 

3.5. Comparison of 15 Mq models 

Comparing Model M15LS-rot with the other 2D sim- 
ulations of our sample allows one to gain insight into 
various aspects that may have important influence on 
the post-bounce evolution. 

3.5.1. High-resolution Model M15LS-rot-hr 

Repeating the computation of Model M15LS-rot dur- 
ing the crucial phase of growing SASI power between 
~420 ms and ^670 ms after bounce with finer radial zon- 
ing (Model M15LS-rot-hr) does not reveal any significant 
differences due to the resolution until 550 ms (cf. Figs. 6- 
9). Only the transient phase of damped SASI activity 
around 500 ms in our reference run of Model M15LS- 
rot (see the middle panel of Fig. 5), which leads to a 
short pause in the expansion of the average shock radius 
(Fig. 6) and a step-like interruption of the rise of the ad- 
vection timescale (Fig. 7), is absent in Model M15LS-rot- 
hr. This suggests that the feature is just a consequence 
of the chaotic and stochastic behavior of the non-radial 
hydrodynamic instabilities that stir up the flow between 
shock front and forming neutron star. 

The conditions that can be considered as favorable 
for an explosion evolve in a promising way also in 
Model M15LS-rot-hr. The average shock radius and the 
mass in the gain layer exhibit a monotonic trend of in- 
crease (Fig. 6 and top panel of Fig. 8, respectively). The 
advection timescale shows the same behavior. Because 
at the same time the neutrino-heating timescale, heating 
rate, and heating efficiency settle to stable levels (Figs. 7 
and 8), the timescale ratio climbs to a value of ~1.3, 
around which it fluctuates at t > 580 ms. Nevertheless, 
after t ~ 570 ms the evolution of Model M15LS-rot-hr 
departs from that of Model M15LS-rot. The latter de- 
velops an explosion after t ~ 570 ms, whereas the former 
reveals SASI activity with a growing amplitude but no 
onset of strong outward shock acceleration until the end 
of our simulation. 

An inspection of the different energies that character- 
ize the energetic conditions in the gain layer reveals that 
the kinetic energy, in particular the part that is associ- 
ated with the radial velocity component, begins to devi- 
ate between the two models aroimd the time when their 
evolution starts to diverge (upper left panel of Fig. 9). 
Also the internal energy of the gain layer begins to grow 
faster in Model M15LS-rot than in its better resolved 
counterpart, but the difference appears slightly later (at 
^^590 ms instead of ~570 ms) and grows more gradually 

^ A longer accretion phase of the nascent neutron star after the 
onset of the supernova explosion is also crucial for setting the ex- 
plosion energy in case of the acoustic mechanism as discussed by 
Burrows et al. (2007c). 
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than for the kinetic energy. Other energies as well as 
the neutrino-heating rate do not show any peculiar dif- 
ferences. While the divergence of the internal energies is 
likely to be a consequence of the stronger shock expan- 
sion in Model M15LS-rot, which leads to growth of the 
mass in the gain layer, the difference of the radial kinetic 
energy could be causal for the discrepant later evolution 
of both simulations. Unfortunately, because of the com- 
plex interdependence of different effects (hydrodynamics, 
gravity, neutrino physics) we are not able to unambigu- 
ously track down the origin of this energy difference and 
to clarify whether it is really the reason or also just the 
result of the evolution that separates at about 570 ms. 

In any case, however, it shows that in the given sit- 
uation the later development of Models M15LS-rot and 
M15LS-rot-hr seems to be extremely sensitive to rela- 
tively small effects. An excess of order 10'*'^ erg in the 
kinetic energy of the radial motions (and also - but a 
factor of ~4 lower - in the lateral kinetic energy, see up- 
per left panel of Fig. 9) and a correspondingly stronger 
SASI and convective activity, which could be caused by 
chaotic fluctuations, appears to be sufficient to initiate 
the outward acceleration of the shock and to cause the 
bifurcation of the subsequent evolution of both models. 
This by itself is a very interesting result, which, of course, 
raises questions about the importance of stochastic vari- 
ations or about the robustness of the explosion found 
in Model M15LS-rot, for example with respect to differ- 
ences in the spatial resolution used in the simulations. 
We repeat and stress here, however, that a variety of di- 
agnostic parameters like the growing shock radius and 
SASI activity as well as the large value of the ratio of 
advection to heating timescale, exhibit an overall trend 
that appears very favorable for an upcoming explosion 
also in the case of Model M15LS-rot-hr. Unfortunately, 
without being able to run this model (and any of the 
other 15 Mq simulations, in particular with high resolu- 
tion) for a sufficiently long post-bounce time — which 
is presently prevented by the limited computer resources 
available to us — we have no convincing support for this 
reasoning. 

3.5.2. Rotating Model M15LS-rot9 

Model M15LS-rot9 was computed with a somewhat 
weaker effective relativistic gravitational potential than 
our standard Model M15LS-rot (Case A of Marek et 
al. 2006 instead of Case R; the former yields the more 
accurate treatment for nonrotating stellar core collapse 
according to comparisons with fully relativistic simu- 
lations). A comparison of these two simulations with 
the nonrotating Model M15LS-2D, which shares the 
description of the effective relativistic potential with 
Model M15LS-rot9 (cf. Table 1), allows us to draw con- 
clusions on the importance of rotation and of the strength 
of the gravitational potential. 

Model M15LS-2D without centrifugal effects but with 
weaker gravity is astonishingly similar to Model M15LS- 
rot in many of the evaluated aspects. This is true, e.g., 
for the neutrino luminosities at t > 200 ms after boimce 
(bottom panel in Fig. 14) and for the neutrino-heating 
rate in the gain layer and the heating efficiency (Fig. 8) 
as well. But also the average shock radii of both models 
are fairly similar and therefore the advection timescales 
do not differ much (Figs. 6, 7). The ratio of the advection 



timescale to the heating timescale shows the same trend 
in both cases, with a slight time lag of Model M15LS- 
2D (lower panel of Fig. 7). This suggests that rota- 
tion does not play an essential role for the evolution of 
Model M15LS-rot towards runaway. In fact, with respect 
to the overall long-time behavior, rotation seems to just 
compensate much of the influence of the stronger gravi- 
tational potential used in Model M15LS-rot. 

The weaker gravity in Model M15LS-rot9 relative to 
Model M15LS-rot allows the shock to move farther out by 
~20 km on average during the first 300 ms after bounce 
(Fig. 6). Since also the infall velocity in the preshock 
region is lower due to the reduced gravitational acceler- 
ation, the advection timescale is significantly longer in 
Model M15LS-rot9. At the same time, the accretion lu- 
minosities of neutrinos in Model M15LS-rot9, in which 
the neutron star is less compact and less hot, are sig- 
nificantly smaller (Fig. 14, bottom panel) and therefore 
the overall postshock heating is weaker (up to roughly a 
factor of two; Fig. 8) and the heating timescale is consid- 
erably longer (typically a factor of two; Fig. 7) than in 
Model M15LS-rot. For this reason. Model M15LS-rot9 
appears to be less favorable for generating a SASI-aided 
neutrino-driven explosion than the nonrotating case with 
the same description of gravity, despite the fact that 
centrifugal support allows for a larger average shock ra- 
dius and correspondingly higher mass in the gain layer 
(Fig. 8). One might therefore conclude that rotation as 
considered in our Models M15LS-rot and M15LS-rot9 
does not improve the conditions for an explosive run- 
away compared to nonrotating models in which low SASI 
modes are able to develop (this comparison turns out to 
be reversed, if such modes are suppressed in the nonro- 
tating case by a constraining choice of the grid setup; see 
Fig. 15 in Buras et al. 2006b). 

A reliable assessment of the influence of rotation on 
the development of an explosion, however, requires a con- 
tinuation of simulations like Model M15LS-rot9 to later 
times. Centrifugal effects were predicted to lead to a re- 
duction of the critical neutrino luminosity for the revival 
of the stalled shock near the polar axis (see the linear 
analysis by Yamasaki & Yamada 2005). The effect grows 
strongly with higher angular momentum of the accreted 
matter. In combination with the fact that rotationally 
deformed neutron stars produce enhanced luminosities 
and/or higher mean energies for neutrinos radiated in the 
polar direction (Janka & Monchmeyer 1989a, b; Kotake 
et al. 2003), this might imply more optimistic conditions 
at very late postbounce times when material with higher 
angular momentum from larger initial radii arrives at the 
accretion shock. 

3.5.3. Nonrotating Models M15LS-2D and M15HW-2D 

The comparison of Models M15LS-2D and M15HW- 
2D yields information about the influence of the neutron 
star EoS, which determines many aspects of the post- 
bounce evolution, most importantly the compactness of 
the nascent neutron star and the properties of the radi- 
ated neutrinos. 

The HW-EoS used in Model M15HW-2D is clearly 
stiffer and the neutron star consequently less compact 
(see Figs. 1 and 6). The neutrino luminosities and 
mean energies are therefore lower (Fig. 1 and the bot- 
tom panel of Fig. 14), which leads to a correspondingly 
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reduced heating rate and heating efficiency in the gain 
layer (Fig. 8), a much longer heating timescale, and thus 
also to a smaller timescale ratio (Fig. 7). 

While in the ID simulation with the HW-EoS the shock 
position follows the behavior of the neutron star radius 
and expands to a larger distance from the center than in 
the ID model with the LS-EoS (see Fig. 1 and Eq. (1)), 
the situation is different for the 2D runs. Here the slosh- 
ing of the shock in Model M15HW-2D reaches a clearly 
smaller amplitude than in Model M15LS-2D during most 
phases of the computed ~400ms of post-bounce evolu- 
tion. This is particularly obvious between 200 and 300 ms 
after bounce and points to a significantly lower strength 
of the SASI activity (for more details, see Marek, Janka, 
& Miiller 2008). An inspection of the kinetic energies in 
Fig. 9 (upper left panel) supports this conjecture, and 
a quantitative evaluation reveals that the I = 1 and 
I — 2 SASI modes reach only 20-30% of the power they 
have in Model M15LS-2D. In fact, the SASI activity in 
Model M15HW-2D grows considerably more slowly and 
the low-mode power is orders of magnitude smaller than 
in Model M15LS-2D until nearly 200 ms after bounce (see 
also Marek et al. 2008). 

This suggests that the compactness of the forming neu- 
tron star has a very strong influence on the growth and 
strength of the SASI and of convective activity in the 
gain layer, either directly by the conversion of gravita- 
tional binding energy to kinetic and internal energy in 
the accretion flow, or indirectly via the effects of larger 
luminosities and mean energies of the neutrinos radiated 
during the period of mass accretion and thus stronger 
neutrino heating. A more rapidly contracting and more 
compact neutron star is obviously more favorable in this 
context. Such a behavior could have different reasons, 
for example a softening of the EoS in the interior of the 
nascent neutron star as a consequence of a microphys- 
ical phase transition that is triggered by evolutionary 
changes (gravitational settling, neutronization, heating 
due to compression and conversion of electron degener- 
acy energy to thermal energy) in the star. We therefore 
emphasize that more favorable conditions for an explo- 
sion during the post-bounce accretion phase of the form- 
ing neutron star neither require a lower value of the in- 
compressibility modulus of the EoS for nearly symmet- 
ric matter at saturation density, nor are they directly 
linked to the EoS properties around bounce, which affect 
the shock formation and prompt shock propagation but 
not necessarily the long-time evolution of the accretion 
shock. 

The gradual growth of relevant parameters (heating 
rate, timescale ratio, kinetic energy in the gain layer, 
etc.) towards the end of the computed evolution ap- 
pears promising for the possibility of an explosion in 
Model M15LS-2D at later times: many properties, of this 
model near the termination of the run, for example the 
neutrino luminosities and neutrino-heating quantities, 
are very similar to those of our explosion model M15LS- 
rot at the same time. Initially Model M15HW-2D looks 
less optimistic than Model M15LS-2D, but there are indi- 
cations that some quantities in Model M15HW-2D might 
reverse their trends in the long run. For example the 
shock radius (Fig. 6), the kinetic energy in the gain layer 
(Fig. 9), the SASI power, and the low- mode shock de- 
formation (Fig. 5) exhibit a strong rise after ~300ms. 



which leads to an increase of the advection timescale 
(Fig. 7) and also of the heating in the gain layer (Fig. 8). 
This improvement could be a consequence of the growing 
compactness of the nascent neutron star. Both simula- 
tions. Model M15HW-2D and Model M15LS-2D, how- 
ever, could not yet be carried on for sufficiently long evo- 
lution periods after bounce to see which differences the 
two employed nuclear EoSs might ultimately could make 
for an explosion at even later times. 

3.6. Neutrino emission 

In Figs. 14 and 15 we present some results on the neu- 
trino emission of our exploding Model M15LS-rot. In 
spite of the deeper effective gravitational potential used 
in this simulation (Case R of Marek et al. 2006), the lu- 
minosities and mean energies of the radiated neutrinos 
at late times after bounce (t > 350 ms) are slightly lower 
than those in Model M15LS-1D. This is again the in- 
fluence of rotation, which compensates some of the con- 
sequences of the stronger gravity in Model M15LS-rot, 
an effect that we have already discussed in the context 
of comparing Models M15LS-rot and M15LS-2D, whose 
neutrino huninosities become very similar after 200 ms of 
post-bounce evolution (see lower panel in Fig. 14). 

The mean energies of the neutrinos radiated from 
Model M15LS-rot (averaged over all observer directions 
relative to the rotation axis), however, do not exhibit the 
crossing of (ep^) and (e,y^) seen in the ID calculations, al- 
though both mean energies get very close (compare the 
middle panel of Fig. 14 with the lower right panel of 
Fig. 1). The reason for this difference can be understood 
from the upper right plot of Fig. 15, which displays the 
mean energies of radiated neutrinos as functions of po- 
lar angle. One can see that in the vicinity of the poles 
(at latitudes 6 < 30° and 6 > 160°) the mean energy 
of muon and tau neutrinos is clearly higher than that of 
electron antineutrinos, while away from the poles both 
are very similar and near the equator their order even 
reverses. This, as well as the trough-like shape of the 
curves (which is also present in the muon and tau neu- 
trino flux densities, see upper left panel of Fig. 15), is 
a consequence of the rotational deformation of the neu- 
trinospheric region as visible in the lower two plots of 
Fig. 15. The superimposed short-wavelength variations 
of the average radiated energies and of the neutrino fluxes 
are caused by local downdrafts of accreted matter near 
the poles and close to the equatorial plane. The neutrino 
flux production associated with such downdrafts is visible 
as bright spots above and around the neutrinospheres in 
the lower left panel of Fig. 15. Near the poles, where the 
downdrafts penetrate deep inward and become very hot, 
neutrinos of all flavors are produced in the tips of these 
flow structures, whereas the lower points of equatorial 
downflows emit only i^e and (this explains the absence 
of a local maximum of the muon and tau neutrino flux 
density around the equator). 

Because of rotation, the density and temperature gra- 
dients at the poles (Fig. 15, bottom right) are much 
steeper and the radial positions of all neutrinospheres 
there are extremely close together^. Most importantly, 

^ In the context of this paper we have adopted Eq. (28) of 
Buras et al. (2006a) for the definition of the energy-averaged neu- 
trinosphere. This means that the average neutrinosphere coincides 
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tlio region where the dominant part of the flux is gener- 
ated before the flux density gets diluted due to scattering 
losses and the growing distance from the source, is sig- 
nificantly below the energy-averaged neutrinosphere for 
muon and tau neutrinos (see the yellow areas in the right 
half of the lower left plot of Fig. 15). This means that be- 
fore the i/f^ and Uj- reach their neutrinosphere, they have 
to diffuse through an extended scattering atmosphere 
that acts like a high-energy filter, i.e., that is inefficient 
in producing more of these neutrinos, but in which en- 
ergy transfers in neutrino-nucleon and neutrino-electron 
scatterings down-grade high-energy neutrinos in energy 
space and reduce the energy flux of heavy-lepton neutri- 
nos (Raffelt 2001; Keil et al. 2003). Moreover, i/g and Pe 
in contrast to heavy-lepton neutrinos are abundantly ab- 
sorbed and emitted by charged-current /3-processes and 
thus are more strongly created in an extended but cooler 
and less dense accretion layer (visible, for example, also 
by a significantly more extended neutrinospheres at low 
latitudes). Both effects in combination are responsible 
for the convergence or even inversion of the mean energies 
of the radiated //,. and in equator-near regions (for a 
discussion of this phenomenon in the context of the non- 
rotating 2D models, see Marek et al. 2008). Close to the 
poles the extreme steepness of the density gradient and 
the additional enhancement of the neutrinospheric fluxes 
of muon and tau neutrinos and antineutrinos by con- 
vection below the neutrinospheres (see next paragraph) 
counteract such a strong filter effect of a scattering atmo- 
sphere. The ID models (cf. Sect. 3.1) define a case that 
is more similar to the equatorial regions of the rotating 
model: although the mass infall rate to the nascent neu- 
tron star drops strongly after the first ~200 ms of massive 
post-bounce accretion, the energysphere of muon and tau 
neutrinos is considerably below the neutrinosphere and 
the density decline below the accretion atmosphere is suf- 
ficiently gradual to establish a thick layer of frequent 
neutrino scatterings between the energysphere and the 
neutrinosphere of heavy-lepton neutrinos. 

Finally, we point out that the region of convective ac- 
tivity inside the nascent neutron star can be discerned 
in the lower left plot of Fig. 15. In the left half of this 
figure, which shows the ffux of z/g, a pattern of green 
and red stripes indicates elongated convective cells that 
are oriented parallel to the rotation axis. This is typical 
of a situation in which convection is constrained by the 
presence of angular momentum gradients (see the discus- 
sion in Buras et al. 2006b). The cell structure becomes 
visible in this image, because Ledoux convection in the 
newly formed neutron star is driven by negative lepton 
number gradients and is therefore associated with vari- 
ations of the neutron-to-proton ratio. Due to the short 
/3-equilibration timescale in the optically thick regime, 
such variations are tightly linked to differences in the 
electron neutrino density as well as flux density. In the 
right half of the lower left image of Fig. 15, one can see 
that a local maximum of the flux coincides with the 
outer edge of the convective region. This suggests that 
convection inside the proto-neutron star causes a notice- 

with the "transport sphere" where the neutrinos in the bulk of the 
spectrum undergo the transition from diffusion to free streaming, 
and not with the "energysphere" where the neutrino production ef- 
fectively ceases and the neutrinos decouple energetically from the 
stellaj: matter. 



able enhancement in particular of the transport of muon 
and tau neutrinos, a fact that was discussed in detail in 
Buras et al. (2006b). 

4. NEUTRON STAR G-MODE OSCILLATIONS 

Recently Burrows et al. (2006, 2007b) discovered in 
their two-dimensional, Newtonian core-collapse simula- 
tions that very late after core bounce the newly formed 
neutron star can be excited to core g-mode oscillations 
with amplitudes of several kilometers. Burrows et al. 
did not obtain neutrino-driven explosions in their sim- 
ulations, but while the shock revealed violent sloshing 
motions due to the SASI instability, the anisotropic gas 
flow around the central object remained gravitationally 
bound and the gas continued to be accreted. After more 
than one second of post-bounce evolution, however, a 
rapid outward acceleration of the shock set in. For an 
11.2M0 progenitor, for example. Burrows et al. (2007b) 
saw this shock revival happening about 1.1 seconds after 
bounce; in a 2OM0 star the outward shock expansion oc- 
curred after 1.2 seconds of post-bounce accretion. In the 
simulations of Buras et al. (2006b), however, the same 
11.2 M0 model was fotmd to explode about 200 ms after 
bounce by the SASI-supported neutrino-heating mecha- 
nism (see Sect. 3.3). The very late explosions in the runs 
of Burrows et al. were powered by the acoustic flux that 
was sent by large-amplitude vibrations of the neutron 
star into the medium surrounding the compact remnant. 
On travelling outward, the pressure waves steepened into 
shocks and dissipated their energy in the medium behind 
the supernova shock, thus creating the thermal energy 
that provided the pressure for the shock acceleration. 

Burrows et al. (2006, 2007b) estimated that the acous- 
tic energy flux generated by this gravity-wave activity 
of the forming compact remnant dominates the energy 
deposition by neutrinos in their models at times later 
than about one second after bounce. A sizable fraction 
of the accretion power is thus converted by the oscillat- 
ing neutron star, which acts as a transducer, to acoustic 
flux instead of neutrino emission. As long as anisotropic 
accretion goes on, the gravity waves are further excited 
and the accretor produces a steady flow of acoustic en- 
ergy. In order to provide a dynamically relevant rate of 
energy input to the shock and possibly the supernova 
explosion, however, the amplitude of the core g-mode os- 
cillations must be sufficiently large at the surface of the 
nascent remnant; several kilometers seem to be necessary 
for creating an acoustic energy flux of a few 10^° erg s"""^ 
(Eq. (1) in Burrows et al. 2007c). 

A closer and independent investigation of this inter- 
esting and potentially important phenomenon is highly 
desirable. Many questions remain to be answered, for 
example: Which decisive factors does the excitation of 
core g-modes due to turbulence and anisotropic accre- 
tion in the SASI layer depend on and how efficient can 
it be? How strong is the acoustic coupling between the 
oscillating neutron star and its surroundings and how 
powerful can the acoustic energy flux become? What is 
the catise of the sudden increase of the g-mode activ- 
ity and of its driving force at late post-bounce times in 
the simulations of Burrows et al. (2006, 2007b)? What 
are the numerical requirements for tracking this phe- 
nomenon in supernova simulations, e.g., is the momen- 
tum conserving treatment of the gravity source term in- 
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troduced by Burrows ot al. important? What is the in- 
fluence of the grid configuration (choice of coordinates) 
and of the grid resolution on the treatment of g-mode ex- 
citation and de-excitation? Does the phenomenon also 
show up in similar strength in 3D instead of 2D sim- 
ulations? In fact, Yoshida, Ohnishi, & Yamada (2007) 
have recently explored the generation of neutron star g- 
modes by external SASI- and turbulence-induced pres- 
sure fluctuations. However, they extracted the boundary 
condition for a semi-analytic g-mode analysis from inde- 
pendent hydrodynamic SASI simulations. While these 
latter simulations were performed with Newtonian grav- 
ity, the g-mode excitation of the neutron star was dis- 
cussed on the basis of a general relativistic stellar model. 
This approach therefore lacks consistency in many re- 
spects. Also, Yoshida et al. (2007) left crucial points un- 
addressed: Docs the prcssiire fluctuation spectrum vary 
with different numerical resolution for the SASI calcula- 
tions? How strongly does neutrino bidk viscosity damp 
the g-mode excitation? Is the assumed boundary con- 
dition appropriate for representing the coupling of the 
neutron star core and the exterior SASI layers, and how 
sensitive is the result to changes of this boundary condi- 
tion? A serious concern of different sort was addressed 
in a paper by Weinberg& Quataert (2008). They calcu- 
lated the damping of the primary I = 1 g-mode in the 
core of the proto-neutron star by its coupling to higher- 
order modes, whose short wavelengths cannot be resolved 
in numerical simulations. They found that the primary 
mode should saturate at an energy (~10^^erg) that is 
much too low to make acoustic power a signiflcant ener- 
getic driver of supernova explosions. 

We do not intend to fundamentally challenge here the 
possibility that the acoustic mechanism might work as 
proposed by Burrows et al. (2006, 2007b, c), because none 
of our simulations was continued to the late times when 
they see the corresponding large-amplitude g-mode os- 
cillations of the neutron star core in their models. In 
the following sections we will concentrate only on two 
questions: Do neutron star g-mode oscillations grow to 
important amplitudes in our simulations and does the 
associated acoustic power radiated by the neutron star 
contribute to the shock revival seen in our models at any 
significant level? 

4.1. Core g-modes in our model runs 

In order to answer these questions, we have performed 
an analysis of the gravity-wave activity in the region 
around and below the neutrinospherc. A direct quan- 
titative comparison of the g-mode strength in the neu- 
tron star core between our simulations and those of Bur- 
rows et al. (2006, 2007b, c) would be desirable and could 
be illuminating, but information suitable for this pur- 
pose is hardly available in the published papers. Fig- 
ure 7 in Burrows et al. (2006), which displays a set of 
time-dependent Z-mode amplitudes of the spherical har- 
monics expansion of the fractional pressure variations 
[P(r, 9) - (P(r, 9))g]/ (P(r, 9))g at a fixed radius of 35 km, 
could be considered as a basis for a straightforward com- 
parison. However, for several reasons the displayed quan- 
tity is not the most appropriate one for diagnosing the 
presence of g-modes in the supernova core and neutron 
star. (1) The mentioned figure provided by Burrows et al. 
was not meant to display a signature of the core g-mode 



activity but was intended to reflect the pressure fluctua- 
tions in the region between neutron star and shock. The 
authors evoked this measure because it is a gauge of the 
potential of pressure fluctuations to excite core g-mode 
oscillations in the more quiescent interior (A. Burrows, 
private communication). The primary origins of these 
"outer" pressure fluctuations are, of course, the SASI and 
neutrino-driven convection, not core pulsations. Without 
a well understood theory of the excitation mechanism 
of the large-amplitude low-order g-modes in the neutron 
star core, however, it is unclear what the external pres- 
sure fluctuations can tell about the fluid motions to be 
expected inside the neutron star. (2) Moreover, pressure 
fluctuations are not an optimal measure of the strength 
of core g-mode pulsations, which rely on buoyancy as 
the restoring force. This is visible, for example, from 
the eigcnfunctions of neutron star core g-modes given 
by Yoshida et al. (2007), where the pressure-dependent 
amplitude appears unspectacular while the Lagrangian 
displacement amplitude can still be large. (3) Pressure 
fluctuations as a function of time at a fixed radius do 
not reflect the evolution of the conditions at a certain 
enclosed mass, but they also include a long-time evolu- 
tionary trend due to the shrinking of the forming neu- 
tron star (and due to its growing centrifugal deforma- 
tion in the case of rotating models). In particular at 
radii close to the boundary layer between the neutron 
star and the violently turbulent, SASI-perturbed post- 
shock region, a gradual increase of the pressure fluc- 
tuations with time does not necessarily signal a rising 
strength of the g-mode driving force at the surface of 
the neutron star core. Instead, the observed growth of 
the pressure fiuctuations can just be a consequence of 
the neutron star contraction, in course of which the ra- 
dial position moves out of the neutron star and deeper 
into the stirred SASI layer. The range of radii that is 
affected by this shift depends on the compactness of the 
neutron star and therefore on the softness of the high- 
density equation of state and on the depth of the gravi- 
tational potential. Since we describe gravitational effects 
by the effective relativistic potential, our neutron stars 
are more compact (in particular with the softer EoS of 
Lattimer & Swesty 1991) than those in the Newtonian 
simulations of the Tucson group. Comparing simulations 
performed with different nuclear equations of state and 
different treatments of gravity (Newtonian or relativis- 
tic) is therefore not conclusive when the same radius of 
evaluation is picked. Simply for this reason a detailed as- 
sessment of the similarities and differences between our 
simulations and those of Burrows et al. (2006) is very 
difficult on grounds of the limited data published. 

Accordingly, an evaluation of the fractional pressure 
variations at different locations in the neutron stars of 
our simulations, guided by Fig. 7 in Burrows et al. (2006), 
leaves too much room for different interpretations and 
does not really illuminate the g-mode activity in the in- 
terior of the nascent neutron star. We therefore decided 
to present a variety of alternative quantities whose be- 
havior can better reflect nonradial fluid perturbations in 
different regions of the supernova core. 

Figures 16 and 17 show the time evolution of low- 
mode spherical harmonics amplitudes (for I = 1,...,5, 
defined in agreement with Eq. (3)) of the velocity per- 
turbations, Vr{r, 6) — {vr{r, 0))g, in the central regions of 



18 



Marek and Janka 



Models M15LS-rot and M15LS-2D, the former figure be- 
tween the neutron star center and a radius of 60 km, the 
latter figure at the positions where the laterally averaged 
density has values of {p)g = 10^^ gcm~^ and 10^^ gcm~^. 
It is obvious that the interior of the neutron star at high 
densities, where the second density location is chosen 
to be, is much more quiet than the environment of the 
neutron star and its surface- near layers, which are rep- 
resented by the first value of the density. The stellar 
plasma in the neutron star surroundings is stirred by the 
violent SASI and convective overturn motions, which in- 
stigate g-modc activity in the outer layers of the compact 
remnant. Figures 16 and 17 reveal that the maximum 
amplitudes of the velocity perturbations deep inside the 
neutron star are of the order of 10^ cms~^, whereas they 
are up to a factor 100 higher in the SASI region and 
adjacent shells. The larger values in Model M15LS-rot 
compared to Model M15LS-2D are a consequence of the 
rotational flattened compact remnant in the former sim- 
ulation, which causes the radius for (p)^ = lO^^gcm"^ 
to reach deeper into the SASI-perturbed layer above the 
poles of the neutron star. 

One should note that the contour for an average den- 
sity of 10^^ g cm~^ is located near the inner edge of the 
convectively unstable shell inside the neutron star, at the 
transition of this layer to the convectively stable region 
closer to the center (roughly at r 10 km). A major 
fraction of the activity seen at this location is therefore 
linked to convection and convective overshooting and not 
to pure gravity waves, but the maximum amplitudes are 
in the range typically found also at positions deeper in- 
side the core. In Model M15LS-2D the dipolc amplitude 
of the velocity perturbations exhibits a characteristically 
different low-frequency modulation, which is absent in 
Model M15LS-rot (see Fig. 17 for {p)g = lO^gcm-^), 
where it is also more difficult to discern the convective 
region inside the neutron star (Fig. 16). This can prob- 
ably be explained by the nature of convection in rapidly 
rotating neutron stars, where due to the Solberg-term 
in the instability criterion the convective cells are con- 
strained to tube-like, narrow structures oriented parallel 
to the rotation axis and located in regions where the an- 
gular momentum gradient is flat (see also Sects. 3.2.1 
and 3.6). The effects of the convection in such an envi- 
ronment are therefore strongly damped compared to the 
nonrotating case. 

Interestingly, at a given density and time, all modes 
from I = 1 to I = a are similarly strong, with a tendency 
of slightly (roughly a factor of two) higher amplitudes for 
Z = 3, 4, 5 (Fig. 17). In the innermost core of the neutron 
star (at r < 10 km) the opposite effect is present and the 
difference of the amplitudes is larger: the amplitudes for 
/ = 1 and / = 2 are the clearly dominant ones during 
most of the post-bounce evolution, I = 3 is intermediate, 
and I = 4, 5 are signiflcantly weaker. 

Consistent with the activity level observed in the veloc- 
ity perturbations, we also find very small fractional den- 
sity variations at given positions, at most some tenths of 
a percent in the neutron star interior and about 10% just 
outside of the neutrinosphere. There is no sign of any sig- 
nificant coherent movement of the neutron star core as 
suggested by Fig. 3 in Burrows et al. (2007c), where the 
isodensity surfaces show large relative displacements of 
their geometrical centers instead of being concentric. 



In order not to just rely on a potentially misleading in- 
terpretation of local fluctuations of the mentioned quan- 
tities, we have also evaluated Models M15LS-rot and 
M15LS-2D for the time-dependent variations of the ra- 
dial positions of certain density values in chosen direc- 
tions (Fig. 18). In the case of a coherent translatory 
movement of some core region, which for symmetry rea- 
sons in the 2D simulations has to occur along the 2;-axis of 
the polar grid, this would also be a rough measure of spa- 
tial shifts of the geometric center of isodensity surfaces. 
Such variations are not found on the scale of the plots in 
Fig. 18 for matter at densities above p > lO^^gcm"'', 
while they are around 0.5 1 km for densities between 
5 X lO^^gcm^^ and lO^^gcm^^ and even larger only 
for lower densities outside of the neutron star^. There, 
however, they originate from the SASI and convective 
overturn activity, which also stirs the neutron star sur- 
face layer below. This interpretation is consistent with 
the fact that one detects radius variations with the same 
behavior and similar or even larger magnitude in equator- 
near regions, where violent downflow activity is also ob- 
served chiring the accretion phase of the stalled shock 
(Sect. 3). A direct inspection of the geometrical centers 
of the isodensity surfaces for the density values consid- 
ered in Fig. 18 reveals full agreement with the message 
and quantitative information in this figure: the geomet- 
rical centers exhibit z-clisplacements relative to the co- 
ordinate center of less than ^0.1 km for densities above 
lO^^gcm""^, at most 0.2km for p = lO^^gcm^"^, some 
0.1km for p = 5 x 10^°gcm~^, and 1-3 km only for 
p < 10^°gcm~^. Taking all together we again con- 
clude that there is no evidence for large-amplitude core 
g-modes in the deep interior of the nascent neutron star 
in our simulations. 

4.1.1. Acoustic energy input to explosion 

Although core g-modcs do not seem to play a signif- 
icant role in the discussed models, the surface gravity 
waves caused by the impact of accretion downflows in 
the outer layers of the neutron star could still be an 
important source of sonic waves and secondary shocks, 
which could contribute energy input to the developing 
explosion. Because of the violently turbulent and time- 
dependent conditions in the neutron star surroundings, 
but also because of the possibility to liberate or absorb 
energy by changes of the composition of the stellar gas, it 
is quite difficult to assess in a quantitatively reliable way 
how important this energy input rate really is. Burrows 
et al. (2006, 2007c) therefore refer to a crude analytic for- 
mula that sensitively depends on the product of several 
factors whose values are not well constrained by the nu- 
merical results (see Eq. (1) in both papers). For checking 
the plausibility of their assumptions they compare their 
estimated numbers with the available rate of energy re- 
lease from gas accretion onto the forming neutron star. 

* It is important to note that in Fig. 18 in contrast to Fig. 17, the 
evaluation is not performed at the position of a laterally averaged 
density value, but at the local density in a chosen direction. In 
the case of a rotationally deformed neutron star the polar radius 
for a certain density is significantly smaller than the radius that 
corresponds to the same value of the lateral density average. The 
former position is therefore close the neutron star surface, while 
the latter position is deeper in the SASI-perturbed surroundings of 
the neutron star. 



Standing accretion shock instability and neutrino-driven explosions 



19 



Our approach here is different. Making use of our hy- 
drodynamical results, wc attempt to estimate the rate of 
acoustic energy transfer to the ejected matter by sound 
waves, Esv,, as the difference of all other energies that 
play a role in the evolution equation for the volume in- 
tegral of the total (i.e., internal plus kinetic plus grav- 
itational) energy in some region behind the supernova 
shock. During the phase of simultaneous accretion and 
outflow of matter (stage (2) described in Sect. 3.4) one 
can write: 



Js\vr>o 



£iiit + P + £kin — 



GMp\ 



C I dVpq^ - (p dnr^Vrpe^^^ 

'Vi JS\Vr>0 



dtEtot + Lo 



(12) 



where for simplicity (and to a very good approximation) 
we assume a spherically symmetric, Newtonian gravita- 
tional potential. In this equation ii^tot is the total energy 
in the volume of integration, Vr the radial velocity com- 
ponent, £int and Ekin the densities of internal and kinetic 
energy, respectively, P the gas pressure, p the mass den- 
sity, (7,y the local rate of neutrino energy deposition, dfir^ 
is the surface element for the surface integrals, and AV 
the volume element for the volume integrals. The surface 
integrals are carried out over a sphere S{r) at radius r, 
which is located behind the shock in all parts. The in- 
tegration includes only material that flows through this 
sphere with positive radial velocity, a constraint that is 
indicated by S\vr > 0. The volume integration, which 
also leads to i?tot, has to be performed over the mat- 
ter that fills the volume between the direction-dependent 
gain radius, Rg{0), and the radius r, and which is ulti- 
mately able to leave the gain layer through the surface 
S(r) with positive radial velocity. The radius r is chosen 
large enough so that the neutrino energy deposition at 
bigger radii is negligibly small. If the volume integra- 
tion was done for all matter in the gain layer, i.e. over 
the whole volume Vg, an extra term would have to be 
introduced into Eq. (12) to account for the energy loss 
associated with matter that is advected inward through 
the gain radius to be accreted onto the neutron star. Be- 
cause of the highly turbulent flow around the gain radius 
and the corrugated structure of the gain surface, it is 
very difficult to evaluate the corresponding surface inte- 
gral. We therefore avoid this term and instead introduce 
the factor ( defining the fraction of the neutrino energy 
that is deposited in the material swept up by the shock 
and ultimately ending up as ejected matter (see Eq. (6) 
and associated discussion) . One can relatively easily ob- 
tain a time-averaged number for ( by comparing the mass 
accretion rate through the shock and the growth rate of 
the mass below the gain radius. 

The first term on the rhs describes the gain of energy in 
the integration volume, the second term the loss of energy 
by cjecta mass Iciaving the surface S{r), the third term 
the total energy input to the ejecta by neutrino heating, 
and the last term accounts for the fact that energy in- 
put by neutrinos (and possibly other sources like sound 
waves) does not only have to yield the positive energy 
carried by the ejecta, but has to provide also the energy 



to overcome the gravitational binding of the progenitor 
gas that enters the gain layer tliroiigli the shock. Refering 
to the detailed arguments given in Sect. 3.4, we calculate 
the internal energy always such that all baryons are as- 
sumed to be assembled in iron-group nuclei. This allows 
us to ignore in Eq. (12) additional terms that describe the 
sizable amounts of energy that can be exchanged between 
the reservoirs of rest-mass and internal energy by pro- 
cesses altering the nuclear composition of the medium, 
i.e. by nuclear burning and photodissociation of nuclei in 
the accretion flow and by nucleon recombination in the 
outflow. The quantity ej^^^^j in the last term on the rhs of 
Eq. (12) is the (negative) specific binding energy (com- 
puted again as the sum of the specific internal, kinetic, 
and gravitational energies) of matter that flows in ac- 
cretion funnels from the shock through the surface S(;r) 
towards the neutron star. This energy with the nu- 
clear binding energy included as mentioned before — is 
typically around —1 to — 2MeV per nucleon during the 
considered accretion phase. It is roughly constant when 
the gas falls inward, because the energy absorbed from 
neutrinos is relatively small when the matter is still lo- 
cated at large distances from the neutron star and is also 
small when the accretion flows move with high veloci- 
ties and thus are exposed to neutrino heating only for a 
very short period of time. Therefore its value is not very 
sensitive to the radius where it is determined". 

We have performed the evaluation for Model M15LS- 
rot in a time interval around 690 ms after bounce near 
the end of our simulation. At this time the mass ac- 
cretion rate through the shock is Mace ~ 0.23 M0S~^ 
and the gas accretion onto the neutron star is still go- 
ing on. A time-averaged fraction of C ^ 0.4-0.5 of 
the infalling mass is fed back into the outflow of ejecta, 
and the conditions in the region between the gain ra- 
dius Rg{0) and the radius r of the sphere S can be as- 
sumed as nearly stationary, in which case dtEtot = is 
a reasonably good approximation. Nevertheless, there 
are still significant time variations in the different terms 
in Eq. (12). For a mass outflow rate through S of 

Mout - C^acc ~ 0.05. ..0.1 MqS-i (- Mgain, the growth 
rate of the mass in the gain layer), we obtain an outflow 
luminosity of Lout (9. ..14) x 10^°ergs~^ and a flow 
rate of gravitational binding energy associated with the 
infalling and reejected mass of L^^^^^ ^ Mout (cbind) 
— (1.5. ..3) X 10'''"ergs^^ (we find an average value of 
(^bind) ^ --1.5MeV per nucleon in the accretion funnels 
of Model M15LS-rot). With a relevant neutrino energy 
deposition rate of CQ^ ^ (10. ..15) x 10^"ergs~-^ the rhs 
of Eq. (12) yields an upper limit for additional energy 
input to the developing explosion by acoustic waves of 
(0.5. ..2) X 10^" ergs^^. The acoustic energy flux originat- 
ing from the violent fluid motions caused by the impact 
of accretion flows near the neutron star surface is there- 
fore at least a factor of ^10 lower than the energy deposi- 
tion by neutrino heating in our simulations. Until 700 ms 
after bounce there is no convincing evidence of any sig- 
niflcant energy transfer to the shock by pressure waves 

^ Bernoulli's theorem implies that the constant energy functional 
includes the specific enthalpy w = (eint + P)/p instead of the 
specific internal energy. Thus the constant quantity is actually 
^bind ~^ ^1 P ^bind' This difference, however, is not essen- 

tial for our discussion. 
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originating from the neutron star. Our 15 Mq model 
M15LS-rot as well as the 11.2 Mq progenitor investigated 
by Buras et al. (2006b) and in this paper (Sect. 3.3) de- 
velop neutrino-powered explosions long before the acous- 
tic mechanism was found to cause an explosion in the 
simulations of Burrows et al. (2006, 2007a, c). 

4.1.2. Momentum conservation 

But is our code actually able to follow the excitation 

and evolution of core g-mode oscillations, in particular 
of Z = 1 type, in which case the gas in the stellar center 
participates in the motion? Since a few radial zones in 
the central < 1.7 km of our numerical grid are treated in 
spherical symmetry (to get around too restricted CFL 
timesteps), one might wonder whether a flow pattern 
with gas motion through the stellar center can be rea- 
sonably well described. Such a spherical core might cre- 
ate severe perturbations, for example because the central 
mass, which cannot move out of the grid center, absorbs 
linear momentum and thus may cause a significant viola- 
tion of momentum conservation. In Fig. 19, left plot, the 
evolution of the total linear momentum of the whole gas 
on the computational grid (about 2Mq, of which about 
three quarters are contained in the neutron star) is shown 
for Model M15LS-rot. The right plot of this figure dis- 
plays the corresponding displacement that the center of 
mass of the gas experiences relative to the grid center in 
the direction of the polar grid axis (which coincides with 
the rotation and symmetry axis of the 2D model) . 

The quantities in Fig. 19 exhibit variability on two 
timescales, which can be easily discerned. A long- 
period modulation with a typical frequency of 30 80 Hz 
(timescales of roughly 15-30 ms) is the consequence of 
the SASI sloshing and convective overturn in the post- 
shock layer (cf. Fig. 5) and the corresponding changes 
of the mass distribution around the neutron star. This 
effect is superimposed by high-frequency fluctuations 
with periods of only milliseconds, which is the typical 
timescale of small-scale mass motions in the close vicin- 
ity of the neutron star and of gravity waves in the neu- 
tron star surface and core. Correspondingly, the linear 
momentum along the z-axis changes its sign with a fre- 
quency of several 100 Hz and the associated net displace- 
ment is tiny. Slightly larger, but still very small is the 
displacement caused by the low- frequency variations: the 
center of mass of the whole gas wobbles around the grid 
center with an amplitude of usually less than 50 meters, 
which is much smaller than the radius of our innermost 
grid cell (300 meters). The separation of the center of 
mass from the grid center grows larger only when the 
outward shock expansion becomes very strong and the 
postshock layer gets inflated with a pronounced global 
deformation at t > 600 ms (see Fig. 4). A clear trend in 
one direction seems to be established towards the end of 
the computational run (with a final displacement, how- 
ever, that is still only insignificantly bigger than the first 
radial grid zone) as a consequence of the clear dipolar 
asymmetry of the accelerating postshock gas, which at- 
tains more monicintum in the southern hemisplier(^ (bot- 
tom right panel in Fig. 4). We observe a very similar 
behavior in the case of our 11.2 explosion model, 
where the z-momentum in the northern hemisphere dom- 
inates when the simulation is terminated. The momen- 
tum asymmetry of the supernova gas could ultimately 



imply a sizable recoil of the neutron star in the oppo- 
site direction, which should drive the neutron star away 
from the coordinate center (possibly leading to numer- 
ical problems), but our simulations are not carried on 
long enough to be conclusive. In Model M15LS-rot the 
velocity of the compact remnant is only ^ 20kms~^ at 
t = 700 ms after bounce, and in the 11.2 A'/q model it is 
6kms~^ at 300 ms after bounce; the final neutron star 
kick will be determined only seconds after the onset of 
the explosion (see Scheck et al. 2006). 

A significant drift of the neutron star away from the 
origin of the grid is likely to lead to numerical problems, 
because the gravity solver as well as the treatment of the 
neutrino transport are written for a centrally condensed 
mass in a polar coordinate grid. Until the end of our 
simulations, however, neither this problem nor violation 
of momentum conservation are a serious concern. 

Another conceivable restriction of our simulations 
might be connected to our use of a relativistic monopole 
potential, whose dominant contribution to the gravita- 
tional potential of the neutron star prohibits an accurate 
response of the gravitational field to the nonradial asym- 
metry of the mass distribution caused by the SASI oscil- 
lations, convective mass motions, and gravity waves. One 
should note here, however, that we still have included the 
higher multipoles of the gravitational potential in their 
Newtonian form, and thus our treatment of the gravi- 
tational effects is certainly much less constraining and 
much less approximative than the so-called Cowling ap- 
proximation. The latter is widely used for discussing 
stellar pulsations in the linear regime, although it ignores 
the influence of the oscillation-induced density variations 
on the gravitational field. 

4.2. Core g-mode tests 

In order to test whether our code can handle / = 1 core 

g-modcs with large amplitudes in the neutron star if they 
were excited (e.g. hydrodynamically by pressure varia- 
tions around the neutron star surface due to anisotropic 
accretion flows), we artificially instigated a large dipole 
g-mode by imposing an / = 1, n = 1 perturbation (i.e., 
we assumed one radial node) of the z-componcnt of the 
velocity field at a time late after bounce (alternatively, 
we also tested an Z = 1, n = 2 mode with two radial 
nodes, obtaining similar results). The velocity field was 
imposed with a constant absolute value on the matter 
within a radius of 50 km. It was chosen such that the 
linear momentum associated with this perturbation was 
zero (the node was therefore initially located at 15 km). 
The once perturbed model was followed in its evolution 
with all (micro)physics being the same as in the long-time 
supernova runs. For exploring the response to g-modes 
with different energies, we varied the amplitude of the 
imposed velocity perturbation by a factor of four, using 
values of 500kms~^ and 2000 km s~^, corresponding to 
kinetic energies of 3.7 x 10^^ erg and 5.9 x 10^® erg, re- 
spectively. 

Figure 20 displays results of two such tests in the 
case of Model M15HW-2D. It shows the time evohi- 
tion of several low-Z mode amplitudes of the spheri- 
cal harmonics decomposition of the radial velocity field, 
{vr{r,6) — {vrir,6))g), at a chosen radius of 10km (left 
panel) and for the I = 1 case in the central region with 
60 km radius (right panel; for comparison we provide 
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tlicrc also the dipolc amplitude of the fractional pres- 
sure variations, (P{r,9) — {P{r,6))g)/ {P)g)- One can see 
that after an initial, short relaxation (caused by the fact 
that our chosen velocity perturbation did not correspond 
to an eigenfunction) an essentially pure dipole mode is 
present with a clean periodicity. The amplitudes of the 
radial velocity variations are a factor of 10-20 larger than 
those that we find for the core g-modes in our supernova 
runs (see Figs. 16 and 17). 

The test model with the lower g-mode energy (which 
is computationally cheaper) was run for many cycles and 
shows a decrease of the oscillation amplitude at a fixed 
radius of 10 km by approximately a factor of two on a 
timescale of about 18ms, but this decay slows down lat- 
eron. In principle, there can be different reasons for this 
damping. The central ID region, though small, might 
brake the motion because there might be friction and 
dissipative losses caused in the 2D flow around the cen- 
tral ID core. This possibility, however, seems to be ruled 
out for the core sizes chosen and the numerical resolution 
employed, because test runs with a core radius of 0.8 km 
instead of 1.7km and with a full 2D treatment (i.e., no 
ID core and thus painfully tiny computational timesteps) 
show the same frequency, amplitude, and damping be- 
havior (Fig. 20, left lower panel). Only when the radius 
of the central ID core is increased to ~3km we can ob- 
serve a slightly faster damping and a gradually decreas- 
ing oscillation frequency. 

Another possible reason for the initial damping is the 
reconfiguration of the mode pattern due to the fact that 
the imposed perturbation does not correspond to an 
eigenfunction. This means that oscillation energy will be 
redistributed within several pulsation periods and there- 
fore amplitude damping is observed in some regions while 
the amplitudes in other regions increase. We find damp- 
ing in the whole volume of the neutron star up to a radius 
of ~30 km, so that the motion outside should be ampli- 
fied. The right panels of Figure 20 indeed reveal such a 
trend. In contrast, hardly any energy is transferred from 
the / = 1 mode to higher modes (in spite of unavoid- 
able numerical coupling), whose amplitudes remain very 
small in the whole perturbed volume of the neutron star 
(Fig. 20, left panel). 

Yet another reason for the observed damping, proba- 
bly enhanced by the just mentioned growth of the g-mode 
amplitude near the neutron star surface, is the acoustic 
energy flux that is sent by the moving neutron star into 
its surrounding medium. Assuming that this ultimately 
is the most efficient channel for energy loss of the ringing 
neutron star, we can obtain an upper limit for the cor- 
responding rate of energy outflow. Since roughly 75% of 
the initial mode energy are lost within 18 ms, we estimate 
an acoustic energy flux of '-~^1.5x 10'"'" erg s^^, comparable 
to the value stated by Burrows et al. (2007c). For the g- 
mode activity and amplitudes observed in our supernova 
runs the energy flux from the oscillating neutron star 
must be expected to be smaller and thus can only make 
a small addition to the neutrino energy deposition be- 
hind the shock, in agreement with our conclusions based 
on the estimates in Sect. 4.1.1. 

The test calculation with the larger excitation am- 
plitude was computationally expensive because of big- 
ger oscillation-induced variations of the thermodynamic 
quantities that led to shifts of the /3-equilibrium between 



neutrinos and the stellar gas. The resulting high neu- 
trino source terms enforced small timesteps in the neu- 
trino transport. The test run was therefore continued 
only for about three full oscillation cycles, but the evolu- 
tion visible in the upper left panel of Fig. 20 agrees nicely 
with the run for the reference case. 

These tests give us confidence that our long-time simu- 
lations should be able to trace also large-amplitude core 
g-modes, if an efficient mechanism was driving their exci- 
tation. We have no reason to suspect that our numerical 
code might fail to describe the hydrodynamic driving due 
to turbulence and anisotropic accretion that Burrows et 
al. (2006, 2007a, c) consider as responsible for instigating 
the large core g-modes in their models (see also Yoshida 
et al. 2007). Therefore we are tempted to conclude that 
until the end of our simulations the conditions for such 
an efficient excitation of gravity waves in the neutron star 
core do not seem to be present in the supernova center. 

5. DISCUSSION AND CONCLUSIONS 

We have presented evidence from our 2D neutrino- 
hydrodynamic stellar core-collapse sinnilations that the 
neutrino-driven mechanism may explain the explosions of 
progenitor stars in a wider range of masses. In this work 
we considered a 15 Mq progenitor (sl5s7b2 of Woosley 
& Weaver 1995) and an 11.2 model (Woosley et al. 
2002). The former is a typical representative of super- 
nova progenitors in the intermediate mass range between 
somewhat more than 10 Mq and around 20 Mq, the lat- 
ter is more typical of stars near the lower mass end of su- 
pernova progenitors with iron cores. A neutrino-driven 
explosion of the 11.2 M© star was found before in a 2D 
simulation by Buras et al. (2006b). Results of a contin- 
uation of this run were also discussed here. Similar to 
the explosion obtained for the 11.2 Mq star, the delayed 
shock revival by neutrino heating in the 15 Mq model 
was fostered and enabled by the presence of a strong, 
low-mode SASI oscillations of the postshock layer^°. The 
SASI modes grow rapidly in the accretion flow to the neu- 
tron star even at conditions where buoyancy instabilities 
are damped in the rapidly infalling gas behind the ac- 
cretion shock and therefore convection in the gain layer 
is initially weak (see Foglizzo et al. 2006, Scheck et al. 
2008). 

The SASI has two crucial consequences, which arc 
favorable for the possibility of neutrino-driven explo- 
sions. On the one hand, the violent sloshing motions of 
the postshock layer with fast expansion and contraction 
phases of the shock create steep entropy gradients in the 
postshock flow, which lead to powerful secondary con- 
vection and thus efficient overturn of the gas in the gain 
region (for a detailed discussion, see Scheck et al. 2008). 
On the other hand, the SASI also pushes the shock to 
a larger average radius. This reduces the mean accre- 
tion velocity in the gain layer (because the velocity in 
the infall region ahead of the shock drops inversely with 
the square root of the shock radius) and thus increases 
the advection timescale of gas from the shock to the gain 

Although our successful explosion for the 15 Mq star was ob- 
tained in a simulation in which we had imposed angular momentum 
on the progenitor model, our comparion with the other simulations 
for this progenitor suggests that the presence of rotation was not 
essential for the development of the explosion (see Sects. 3.5.2 and 
3.5.3). 
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radius. The accreted matter is therefore exposed to neu- 
trino heating for a longer time, a fact that is reflected by 
a significant increase of the critical ratio of the advection 
timescale to the neutrino-heating timescale. A strong 
SASI activity has the consequence that this timescale 
ratio comes much closer to the value of unity, signaling 
favorable conditions for a neutrino-driven explosion. The 
combination of these effects indeed turned out to drive 
the 15 Mq model (specifically our Model M15LS-rot) to- 
wards a runaway instability at about 600 ms after core 
bounce, while for the 11.2 Mq star this happened roughly 
200 ms after bounce. 

Unfortunately, due to the considerable demand of com- 
puter time with the high numerical resolution needed 
for converged results, we were so far not able to con- 
tinue our simulations to the phase where the blast is 
fully developed and the explosion energy can be deter- 
mined from our computations. However, several factors 
are indicators that the 15 Mq Model M15LS-rot consid- 
ered here just like the investigated 11.2 model is in 
the process of undergoing the transition from collapse 
to successful outburst. This is suggested not only by 
the rapidly accelerating expansion of the shock radius 
near the end of the simulations and by the fact that the 
critical timescale ratio exceeds unity permanently with 
rapidly growing value. It is also suggested by a continu- 
ous trend of increasing energy in the gain layer. Before 
Model M15LS-rot was stopped, a significant and growing 
gas mass in the gain layer had obtained a positive spe- 
cific energy and was thus ready to become gravitationally 
unbound and to promote further outward propagation of 
the shock against the gravitational pull of the neutron 
star. The mean shock radius has already reached about 
600 km. This is sufficiently large for the temperature in 
the postshock ejecta to be so low that considerable nu- 
cleon recombination to a-particles has set in and for the 
gas swept up by the shock to be so cool that heavy nu- 
clei do not experience complete dissociation any longer. 
The release of nuclear binding energy and the reduction 
of dissociation energy allow for a higher pressure behind 
the shock and support the shock expansion. Last but 
not least, the shock has entered the progenitor layer (at 
about 500 km) where silicon is not yet burned to iron, 
and its most extended parts reach the shell where oxygen 
is still unburned (at ~700km). According to Bruenn et 
al. (2006) and Mezzacappa et al. (2007), shock- initiated 
nuclear burning and the reduced ram pressure in these 
layers will assist the ongoing strong neutrino heating in 
driving the explosion. For all these facts we do not sec 
any reason to suspect that the outward acceleration of 
the shock near the end of our simulated Model M15LS- 
rot might break down again. 

The energy of such SASI-supported, neutrino-driven 
explosions is not determined at the moment when the 
runaway sets in. Instead, accretion of gas towards the 
gain radius and the neutron star can proceed at the same 
time as shock acceleration sets in due to the expansion of 
rising, neutrino-heated, buoyant gas. This characterizes 
the generically multi-dimensional nature of the develop- 
ing explosion. The ongoing accretion continuously chan- 
nels fresh material, which is swept up by the shock, to 
the region near the gain radius, where neutrino heating 
is strongest. Roughly half of the gas penetrates to the 
cooling layer below the gain radius and ultimately settles 



onto the neutron star, the other half absorbs with high 
efficiency energy from the neutrinos radiated by the neu- 
tron star. This gas turns around and rises again, driven 
by buoyancy forces, and contributes to the energy of the 
explosion. Our estimates show (and hydrodynamic sim- 
ulations by Scheck et al. 2006 have demonstated) that 
simultaneous accretion and shock acceleration can go on 
for many hundred milliseconds and therefore the final ex- 
plosion energy can only be determined by computations 
that cover a much longer time evolution than we were 
able to follow in our current simulations. Our analytic 
considerations, based on the situation at the beginning of 
the explosion, reveal that the 11.2 M© model can be ex- 
pected to develop an explosion energy of several 10^° erg, 
while due to the higher mass accretion rate the canonical 
value of 10^^ erg is well in reach of the considered 15 Mq 
star. The final explosion energy should increase sensi- 
tively when the initial expansion of the shock is relatively 
slow and when the progenitor maintains an appreciable 
mass accretion rate for a longer time. 

Because the onset of the discussed SASI-supported 
neutrino-driven explosions for 11.2 and 15 Mq stars is 
considerably delayed, the neutrino-heated gas mass in 
the gain layer at the time when the blast takes off is 
fairly low and therefore a powerful gas ejection has to 
rely on a longer lasting phase of simultaneous accretion 
and outfiow. The possibility of such stable accretion with 
downflows from one (or more) direction(s) and gas out- 
flow in the other direction(s) was indeed seen in the sim- 
ulations of Scheck et al. (2006). The significance of this 
generically multi-dimensional phenomenon in powering 
the explosion to completion and in achieving the neces- 
sary energies contradicts the concept of supernovae be- 
ing energized by a spherically symmetric neutrino-driven 
wind as suggested by Burrows & Goshy (1993). This has 
been realized recently also by Burrows et al. and was ad- 
vocated by them in the context of the acoustic explosion 
mechanism (Burrows et al. 2007c), which requires the 
accretion of gravitational/mechanical energy, and in the 
context of magnetohydrodynamical explosions (Burrows 
et al. 2007a), which require the accretion of differential 
kinetic energy. 

The strong shock acceleration, which we found to be- 
gin about 200 ms after bounce in the 11.2 Mq progenitor 
and at 600 ms after bounce in the investigated 15 M© 
model, is associated with a large dipolar deformation of 
the shock. This is a clear indication for the dominance of 
the lowest nonradial SASI mode at the onset of the run- 
away. It suggests that the developing explosion is likely 
to be very asymmetric, similar to the artificially initi- 
ated neutrino-driven explosions that were computed by 
Scheck et al. (2006) and led to neutron star kick veloc- 
ities (with maximum values of more than 1000 km s^^) 
in agreement with measurements. Also the large-scale 
anisotropics and mixing processes that can be found in 
many supernovae and that seem to determine the appear- 
ance of their gaseous remnants might be a consequence of 
the pronounced asymmetry of the developing blast waves 
in our simiflations (see Kifonidis et al. (2006) for studies 
of the long-time evolution of anisotropic neutrino-driven 
explosions). 

We compared 2D models with and without rotation 
and with a stiff and a soft nuclear equation of state More- 
over, we tested the influence of different prescriptions of 
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the effective relativistic gravitational potential in our ap- 
proximation to fully relativistic calculations. We found 
that a nascent neutron star that contracts faster and 
is more compact, which typically happens in the case 
of a softer high-density equation of state (or a stronger 
gravitational potential) , has a favorable influence on the 
SASI-supported neutrino-driven mechanism. This is a 
consequence on the one hand of a more rapid growth 
of low SASI modes, amplified by the release of gravita- 
tional binding energy, and on the other hand of higher 
luminosities and mean energies of the emitted neutrinos 
during the accretion phase and correspondingly stronger 
neutrino heating in the gain layer. The latter also sup- 
ports the development of violent SASI and convective 
activity in the postshock region. We stress that the EoS 
properties that lead to a faster compactification of the 
nascent neiitron star during its post-bounce accretion do 
not need to be linked to the EoS characteristics at core 
bounce or to the value of the incompressibility modulus 
of symmetric nuclear matter at saturation density. The 
evolutionary changes in the interior of the remnant could 
also trigger a softening of the EoS by a reduction of the 
adiabatic index, e.g., as a result of a phase transition. 

Rotation, in contrast, can be diagnosed to have the 
opposite effect. Although centrifugal support leads to a 
larger average shock radius, a longer advection timescale 
of the matter falling through the gain layer, and a re- 
duced binding energy of accreted matter, these helpful 
effects are more than compensated by lower luminosities 
and mean energies of the neutrinos radiated by the signif- 
icantly more extended and cooler rotating proto-neutron 
star. The overall neutrino heating and heating efficiency 
found in rotating models is considerably lower and the 
neutrino heating timescale in the gain layer correspond- 
ingly longer. The critical ratio of advection timescale 
to heating timescale therefore signals less favorable con- 
ditions for an explosion in comparison to a nonrotating 
model with otherwise the same physics. This pessimistic 
judgement of the role of rotation is opposed by the as- 
sessment of Yamasaki & Yamada (2005), who found by 
linear analysis that sufficiently rapid rotation leads to a 
reduced critical neutrino luminosity for shock revival in 
the polar direction. Simulations for longer postbounce 
times are needed to obtain conclusive information about 
this possibility. 

Although we clearly see the presence of gravity-waves 
in the newly formed neutron star, we cannot detect any 
important influence on the dynamical evolution of the be- 
ginning supernova explosion. On the one hand, the am- 
plitudes of core g-mode oscillations in our supernova runs 
remain very small. On the other hand, the estimated 
flux of acoustic power associated mainly with the con- 
siderable surface gravity-wave activity is dwarfed by the 
rate of the neutrino energy deposition. The small ampli- 
tudes of core g-modes found in our models are not in dis- 
agreement with the calculations of Burrows et al. (2006, 
2007b, c), who discovered violent neutron star core mo- 
tions only at very late times {t>ls) after core bounce, 
which we are unable to reach in our simulations. 

We convinced ourselves by numerical tests with artifi- 
cially instigated large-amplitude core g-mode vibrations 
that our supernova code is well able to capture this effect 
even in the case of sizable core displacements (of order 
one kilometer and more). Moreover, we do not see any 



reason to suspect that our numerics could be unable to 
track the g-mode excitation by pressure fluctuations due 
to anisotropic accretion and turbulence in the SASI layer, 
in particular since we make sure to have good numerical 
resolution inside the neutron star as well as in the re- 
gion between the neutron star and the shock. Both facts 
together make us confldent that we should observe large- 
amplitude core pulsations if the physical conditions were 
present to drive their excitation. 

The 11.2 M0 and 15 Mq simulations discussed in the 
present work suggest that — at least in this mass range of 
progenitor stars — the SASI-aided neutrino-driven mech- 
anism can lead to explosions significantly earlier than the 
acoustic mechanism proposed by Burrows et al. (2006, 
2007b, c), a possibility that is not challenged by the lat- 
ter authors. Naturally, our simulations, which had to be 
terminated at latest ~700 ms after boimce, do not allow 
us to make any statement about the neutron star g-mode 
activity at even later times (with the potential to pump 
additional energy into the already launched blast), nor 
can we exclude a potentially important role of the vibrat- 
ing neutron star as an acoustic energy source for trig- 
gering the explosion in cases where the neutrino-driven 
mechanism is too weak and the explosion is as delayed 
as observed by Burrows et al. (2006, 2007b, c). 

Neutrino- driven explosions have recently been also 
found in 2D simulations with sophisticated, energy- 
dependent neutrino transport by Bruenn et al. (2006) 
and Mezzacappa et al. (2007). Their calculations for 
IIM0 and 15 Mq stars, which in contrast to ours were 
Newtonian and were performed with a flux-limited neu- 
trino diffusion scheme, revealed the initiation of an ex- 
plosion at the time when the inner edge of the oxygen 
layer accretes through the shock. The authors reported 
that explosions were only obtained when they used an al- 
pha network of nuclei for following composition changes 
by nuclear burning in the collapsing stellar layers, but 
not when they applied a "flashing" treatment with in- 
stantaneous conversion of the nuclear abundances from 
non-NSE to NSE. They diagnosed that effective energy 
release by oxygen burning in the immediate vicinity of 
the shock, which in their case happens only with the al- 
pha network, assists the explosion, in particular in case 
of a weak shock. Their papers do not provide the infor- 
mation that would allow us to conduct a detailed com- 
parison of the composition evolution in their simulations 
and in ours. However, we emphasise that although we do 
not use a nuclear reaction network, we still follow compo- 
sition changes in the non-NSE regime by converting oxy- 
gen to silicon and silicon to iron when the corresponding 
burning temperatures are reached due to compressional 
(or shock) heating in the infalling stellar layers. In our 
calculations for the 11.2 Mq and 15 Mq progenitors a va- 
riety of important indicators signal the onset of the blast 
and strong outward shock acceleration already long be- 
fore the accretion shock has reached the inner bound- 
ary of the oxygen-rich layer. This, of course, does not 
mean that oxygen burning in the collapsing layers has no 
bearing on the possibility of an explosion, but it means 
that the inauguration of the blast is not correlated with 
oxygen combustion happening at the shock. Indepen- 
dent of whether a more accurate network treatment of 
nuclear reactions and of the transition to NSE instead 
of our simplified description can make a difference for 
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the shock propagation or not, the simulations by Bruenn 
et aL (2006) and Mezzacappa et al. (2007) suggest that 
any such improvement is hkely not to disfavor explosions 
of the investigated stars. In this sense the fundamental 
agreement of the outcome of these different simulations, 
in spite of various differences in numerical aspects and 
input physics, may be considered as very encouraging. 

Certainly, our current simulations can only be sugges- 
tive for how neutrino heating and hydrodynamic insta- 
bilities in collapsing stellar cores can collaborate to ini- 
tiate the supernova explosion. The simulations must be 
continued to later times for determining the explosion 
properties, and definitely more calculations are needed 
to obtain a clear understanding how the shock revival de- 
pends on the core structure of the progenitor stars and on 
the physics that plays a role in the supernova core, e.g., 
the neutrino opacities, equation of state, rotation, the 
depth of the gravitational potential, and also stochastic 
elements like chaotic fluctuations of the SASI strength in 
the nonlinear regime. Also the influence of the numerical 
resolution has to be more closely studied. Our simula- 
tions cannot assess all the interesting questions involved. 
The results as presented here should only be seen as an 
indication that the neutrino-driven mechanism is a viable 
possibility for driving supernova explosions of progenitor 
stars significantly more massive than lOM©. 

The positive trend towards a SASI-supportcd, 
neutrino-driven runaway late after bounce might be 
linked to a (relatively) soft equation of state for neutron 
star matter, like the EoS of Lattimer & Swesty (1991) 
used by us. Our set of calculations for the considered 
15 Mq star, although including also a model with a sig- 
nificantly stiffcr nuclear EoS, does not provide an an- 
swer to the question how robust the success of the SASI- 
supported neutrino-heating mechanism is to variations of 
such an important physics ingredient. Systematic long- 
time simulations for different nuclear equations of state 
are needed to explore the differences at times later than 
covered by our present calculations. Also the exact mo- 
ment when the explosion sets in can well be sensitive to 
various uncertain aspects like the amount of rotation in 
the stellar core, the possible influence of amplifying mag- 
netic fields, imsettled details of the subnuclcar equation 
of state and of neutrino-matter interactions, neutrino os- 
cillations, the missing treatment of lateral neutrino fluxes 
in our simulations, and the approximative description of 
the effects of rclativistic gravity. Even stochastical vari- 
ations at some fairly low energy level (C(10^'^) erg in the 
gain layer) seem to be sufficient to foster an explosion at 
an earlier time. Simulations in full relativity are certainly 
desirable, in particular for stellar core collapse with ro- 
tation, because in this case the optimal choice of the ef- 
fective gravitational potential is not obvious (see Marek 
et al. 2006). Moreover, supernova modeling in three spa- 
tial dimensions must be the ultimate goal, because the 



growth rate and properties of the hydrodynamic insta- 
bilities in 3D must be expected to differ from the axially 
symmetric 2D case, and new degrees of freedom may play 
a non-negligible and potentially helpful role, for exam- 
ple spiral waves {m ^ modes) and triaxial instabilities 
(see, e.g., Blondin k Mezzacappa 2007, Ott et al. 2007, 
Yamasaki & Foglizzo 2007, Iwakami et al. 2008). 

In the light of our present results, supernova explo- 
sions appear to be an accretion instability rather than an 
aspherical wind as described by Burrows et al. (2007c). 
The SASI and convective instabilities, and in particu- 
lar the simultaneous presence of accretion and outflow 
of neutrino-heated gas, which drives the shock expan- 
sion, are crucial ingredients of the explosions described 
in this paper and constitute the explosions as a gener- 
ically multi-dimensional phenomenon. The gas outflow 
during the main phase of shock revival and the build- 
up of the explosion energy are fed by neutrino-heated, 
accreted gas. Only much later, after the accretion has 
ceased, does the gain radius retreat to the neutron star 
surface and the cooling layer between neutrinosphere and 
gain radius shrinks to a very narrow region or disappears 
completely. This marks the onset of the neutrino-driven 
baryonic wind phase, in which the dilute gas outflow from 
the nascent neutron star is determined solely by the con- 
ditions at the surface of the hot, compact remnant and 
not by the (accretion) properties of the dying star. 
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Table 1. Investigated 15 Mq models. 
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''Number of radial grid points, increasing with time. 
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Number of energy grid points for the neutrino transport. 
''Effective relativistic potential; "old" corresponds to definition introduced by Rampp & Janka (2002), which is equivalent to Case R of 
Marek et al. (2006); "new" corresponds to Case A of Marek et al. (2006). 
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Fig. 1. — Results of our two ID models, M15LS-1D (bold lines) and M15HW-1D. Time is normalized to the moment of core bounec. 
Top left: Central densities as functions of time. The horizontal dash-triple-dotted line marks the value of the nuclear saturation density, 
po ^ 2.7 X lO^'^gcm"^. Top right: Radii of shock (solid lines) and electron neutrinosphere (dashed lines) as functions of time. The 
transient shock expansion around 170 ms after bounce is caused by a composition interface with sudden decrease of the mass accretion 
rate arriving at the shock. The mass accretion rate is displayed by the bold dotted line. Bottom left: Luminosities of electron neutrinos, 
electron antineutrinos, and muon or tau neutrinos (or their antiparticles) as functions of time. Bottom right: Mean energies of radiated 
neutrinos (computed as the ratio of energy to number flux) for electron neutrinos, electron antineutrinos, and heavy-lepton neutrinos and 
antineutrinos. The luminosities and mean energies axe shown as measured by an observer at rest relative to the stellar center at 400 km 
(from where the gravitational redshift to infinity is negligibly small). The softer LS-EoS leads to a much larger central density, a significantly 
more compact neutron star and therefore considerably higher luminosities and mean energies of the radiated neutrinos. 
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Fig. 2. — Time evolution of Model M15LS-rot visualized by mass shell trajectories. In this 2D simulation with rotation, the mass-shell 
lines mark the radii of spheres that contain certain values of the rest mass (the plot is based on an evaluation of the mass-weighted lateral 
average of the 2D data set). They are spaced in steps of 0.025 Mq with bold lines every 0.1 Mq. The thick solid line starting at t = 
denotes the mass-averaged shock position, the blue lines represent the mean neutrinospheres of i/e (solid), i>e (dashed), and heavy-lepton 
neutrinos (dash-dotted), the black dashed curve shows the mean gain radius, and the location of the composition interface between the 
silicon shell and the oxygen-enriched Si-layer of the progenitor star at 1.42 Mq is highlighted by a red dashed line. Different shadings 
indicate regions with different chemical composition. Dark grey marks the layer where the mass fraction of oxygen is larger than 10% 
(which corresponds to the inner boundary of the layers that contain significant amounts of oxygen), medium grey the region where the 
mass fraction of heavy nuclei with mass numbers A > 56 exceeds 70%, the yellow band in between is the layer where both abundance 
constraints are not fulfilled (in this region silicon and sulfur are abundant), light grey indicates those regions where more than 30% of 
the mass is in a-particles, and the white areas enclosed by the shock front contain mostly free nucleons and only a small mass fraction 
(less than 30%) of o-particles. At times t > 600 ms post bounce, slightly darker grey patches in the light-grey postshock regions contain 
a mass fraction of more than 60% helium nuclei. This signals that the nucleon recombination becomes more complete and/or that the 
dissociation of alpha particles to free nucleons is less complete in the matter expanding behind the outgoing shock because of low postshock 
temperatures when the shock reaches larger radii. Note that compressional heating triggers nuclear burning (described in our simulations 
by a "flashing treatment", see Sect. 2.1) and leads to changes of the chemical composition in the infalling stellar layers. 
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Fig. 3. — Radial positions of the shock near the north and south poles of Model M15LS-rot as functions of post-bounce time (white 
lines). The color coding represents the entropy per nucleon of the stellar gas. The quasi-periodic, bi-polar shock expansion and contraction 
due to the SASI can be clearly seen. 
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Fig. 4. — Six snapshots from the post-bounce evolution of Model M15LS-rot. The color coding represents the entropy of the stellar gas. 
The shock is visible as deformed sharp discontinuity between low-entropy, infalling matter in the upstream region and high-entropy, boiling 
matter behind the shock; its position is highlighted by a bold, solid white contour. The top left plot shows the entropy distribution at 
t = 119 ms after bounce, about 40 ms after the postshock convection has reached the nonlinear regime and the shock develops first small 
nonsphericities. The top right and middle left plots {t = 454 ms and and 524 ms after bounce, respectively) demonstrate the presence of very 
strong bipolar oscillations due to the SASI, the middle right plot {t = 610 ms p.b.) displays the beginning of a rapid outward expansion, 
and the lower two plots (for t = 650 ms and 700 ms post bounce) show the onset of the explosion with a largely aspherical shock that 
possesses a dominant Z = 1 deformation mode. Note that the radial scale was adjusted in the last three snapshots and that the contracting 
nascent neutron star exhibits a growing prolate deformation because of the rotation considered in this simulation. The thin, solid white line 
in each panel marks the direction-dependent location of the gain radius, and the thin dotted, dashed, and dash-dotted white lines indicate 
the inner boundaries of the regions where iron-group elements, silicon, or oxygen, respectively, dominate the composition (the contours 
are defined by mass fractions of 30% iron-group elements, 30% silicon, and 10% oxygen, respectively). In some of the panels not all these 
composition interfaces are located within the plotted area, and the iron-dissociation line or the iron-silicon interface can (at least partly) 
overlap with the shock contour. We point out that the rotation of the model is so slow that the composition interfaces in the preshock 
region exhibit no visible centrifugal deformation. 
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Fig. 5. — SASI and convective activity in Model M15LS-rot versus post-bounce time. Top: Coefficients of tlie dipole (Z = 1) and 
quadrupole (Z = 2) modes of the splierical harmonics expansion of the angle-dependent shock position, normalized to the amplitude of the 
I = mode. For better visibility the curve of the 1 = 2 mode is shifted downward by 0.7 units. The results for Models M15LS-2D and 
M15HW-2D are also shown for comparison. Middle: Power of the I = 1 and I = 2 modes of the spherical harmonics expansion of the 
fractional pressure variations [P{r,9) — (f (P)g^ integrated over the volume between average electron neutrinosphere and average 

shock radius. Bottom: Frequency spectrum versus time of the combined power in Z = 1 and Z = 2 modes shown in the middle panel. 
The fourier analysis was performed every two milliseconds, sampling information in time windows of 64 ms width (which leads to discrete 
frequencies of 15.6 Hz and multiples), and the spectra were normalized on each time slice. 
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Fig. 6. — Average shock radii (solid) and average neutron star radii (dashed) for the set of computed 2D models with the 15 Mq 
progenitor as functions of post-bounce time (left vertical scale of the figure). The shock position is defined as the arithmetical average over 
all directions, and the mean radius of the (in case of rotation, deformed) neutron star is determined as the arithmetical mean of all radial 
positions where the density equals 10^^ gem"''. The bold black curve shows the time-dependent mass accretion rate just ahead of the shock 
in Model M15LS-rot (scale on the right vertical axis of the plot). The additional black line between 420 and about 670 ms corresponds 
to Model M15LS-rot-hr, which is a test calculation of Model M15LS-rot with significantly increased radial resolution in the neutron star 
surface and neutrino-heating layer. 
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Fig. 7. — Dynamically important timescales for the set of computed 2D models with the 15 Mq progenitor as functions of post-bounce 
time. Top: Mean timescale for advection of accreted matter from the shock to the gain radius. Middle: Neutrino-heating timescale of 
matter in the gain layer. Bottom: Ratio of advection to heating timescale. The additional black line between 420 and about 670 ms 
corresponds to Model M15LS-rot-hr, which is a test calculation of Model M15LS-rot with significantly increased radial resolution in the 
neutron star surface and neutrino-heating layer. Note that towards the end of our explosion simulation (Model M15LS-rot) the advection 
timescale and the timescale ratio increase steeply because a dominant fraction of the mass in the gain layer begins to expand outward 
behind the shock instead of falling towards the gain radius (this is also reflected by a steep growth of the mass in the gain layer, see Fig. 8). 
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Fig. 8. — Neutrino heating conditions in the gain layer for the set of computed 2D models with the 15 Mq progenitor. Mass {top), total 
net rate of neutrino heating (middle), and neutrino heating efficiency (bottom) as functions of post-bounce time. The net rate of neutrino 
heating is defined as the difference of the neutrino energy deposition rate in the gain layer and the energy loss rate by the reemission of 
neutrinos. The heating efficiency is computed as the ratio of the total net heating rate to the sum of Ue and Pe luminosities as seen in the 
observer frame. The additional black line between 420 and about 670 ms corresponds to Model M15LS-rot-hr, which is a test calculation 
of Model M15LS-rot with significantly increased radial resolution in the neutron star surface and neutrino-heating layer. 
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Fig. 9. — Evolution of energies in the 15 Mq models. The upper plots display integrated energies in the gain layer as functions of 



post-bounce time for several models of our set of 2D simulations, 



the lower plots show energies per nucleon and corresponding 



in the gain layer as functions of time for Model M15LS-rot. Top left: Kinetic energies E'^ 
radial and in radial plus lateral direction, respectively. Top right: Internal energy and rotational energy. Bottom left: Mean total (kinetic, 
including rotational, plus internal plus gravitational) energy per baryon (thick line) and energy range that contains 90% of the mass in 
the gain layer of Model M15LS-rot (dotted). Bottom right: Masses in the gain layer of Model M15LS-rot with total energies per nucleon 
above certain values (left panel) and increase of the "explosion energy" of the model, defined as total energy of all matter in the gain layer 
with positive radial velocity (right panel). In the upper two plots, the additional black line between 420 and about 670 ms corresponds 
to Model M15LS-rot-hr, which is a test calculation of Model M15LS-rot with significantly increased radial resolution in the neutron star 
surface and neutrino-heating layer. The rise of the curves for Model M15LS-rot towards the end of the simulation is a signal of the beginning 
explosion, which is also accompanied by a steep increase of the mass in the gain layer (see Fig. 8). 
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Fig. 10. — Same as Fig. 2 but for our 2D explosion simulation of an 11.2 Mq progenitor star. Note that the mass-shell spacing outside of 
the red dashed line at an enclosed mass of 1.25 Mq (marking the composition interface between the silicon layer and the oxygen-enriched 
Si-shell) is reduced to steps of 0.0125 Mq instead of 0.025 Mq. 
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Fig. 11. — Left panel: Mean shock radius (arithmetical average over all lateral directions, dashed line) and maximum and minimum shock 
positions as functions of post-bounce time for our 2D explosion simulation of an 11.2 Mq progenitor. Right panel: "Explosion energy" of 
the 11.2 A/q star, defined as the total energy (internal plus kinetic plus gravitational) of all mass in the gain layer with positive radial 
velocity, as function of post-bounce time. 
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Fig. 13. — Mass {top left), neutrino-heating rate [top right), heating efficiency {bottom left), and heating and advection timescales {bottom 
right) in the gain layer as functions of time for our 11.2 Mq explosion model. 
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Fig. 14. — Luminosities (top) and mean energies (middle) of radiated neutrinos in Model M15LS-rot, and total luminosities (summed for 
neutrinos and antineutrinos of all flavors) for all 2D simulations with the 15 Mq progenitor (bottom) as functions of time after bounce. All 
quantities are evaluated at a radius of 400 km for an observer in the rest frame of the stellar center, and the average energies are defined 
as ratio of energy flux to number flux. The decline of the curves near the end of the simulation of Model M15LS-rot signals the onset of 
the explosion, which loads to reduced accretion luminosities of Ue and Pe from the nascent neutron star. 
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Fig. 15. — Top: Neutrino energy flux densities {left) and mean energies {right) as functions of polar angle at 540 ms after bounce for 
Model M15LS-rot. The quantities are measured at 400 km and given in the observer frame, and the average energies are defined as ratio 
of energy flux density to number flux density. Bottom, left: Energy flux densities of electron neutrinos and (one kind of) heavy-lepton 
neutrinos in and around the nascent neutron star. The mean neutrinospheres, which are marked by white lines for (solid), Ve (dotted), 
and (dashed), are significantly deformed because of the presence of rotation. Bottom, right: Density and temperature in the vicinity of 
the forming neutron star. 
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Fig. 16. — Amplitude of the I = 1-mode component of the spherical harmonics decomposition of the velocity variations normalized by 
lO^cms"-"^, {vr(r,9) — {vr(r, 9))g)/10^' cm/s, as a function of radius and post-bounce time in our simulation of the rotating 15Mq model 
M15LS-rot (left) and in the corresponding simulation without rotation, M15LS-2D (right), once for the whole post-bounce evolution {top) 
and another time for the last ~100ms of both simulations {bottom). In both cases the central region with a radius of 60 km is shown, but 
the color scale for the amplitude is constrained to a different range of values (much more narrow than the absolute maxima and minima) 
in order to visualize activity in different regions. The interior of the neutron star at r < 10 km is much more quiet than the outer layers 
where SASI and convective motions perturb the surroundings of the proto-neutron star and stir g-mode activity in its surface layers. In 
particular, there is no sign of any sizable core g-mode oscillations in the neutron star core. The dashed lines mark the positions where the 
laterally averag ed density has values of {p)g = 10^'^, 10^^, 10^^, lO", and IQlOgcm-^ (from bottom to top). 
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Fig. 17. — Amplitudes of the I = 1, 2, 3, 4, 5-modes of the spherical harmonics decomposition of the velocity variations, {vr{r,6) — 
{vrir,0))g), as functions of post-bounce time at fixed values of the laterally averaged density (see the dashed lines in Fig. 16) for our 
Models M15LS-rot (left) and M15LS-2D (right). The upper panels show the amplitudes at a density of {p)g = lO^^gcm"^ with a 
normalization by lO^cms"^ (left) and 10*cms~^ (right), the lower panels at {p)g = lO^^gcm"^ with a normalization by lO^cms"^. 
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Fig. 18. — Radius variations rp(t) as functions of post-bounce time of the locations for different values of the local density (as specified in 
the plots) at a position near the pole of the spherical coordinate grid {upper panels) and at a latitudinal angle of 45 degrees {lower panels) 
in the nascent neutron star of Model M15LS-rot {left) and Model M15LS-2D (right). The displayed variations are indicative for major mass 
motions and mass displacements associated for example with convective overturn or g-mode oscillations. One can see that large amplitudes 
are present only in the outermost layers of the neutron star at densities below some 10^" gem"'' (see also the zoom in the inset), where the 
violent SASI and convective activity in the postshock layers around the neutron star makes an impact. The interior of the neutron star is 
essentially quiet, in particular there is no sign of any sizable core g-mode oscillations or pulsational motions of the neutron star core. 
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Fig. 19. — Momentum conservation in Model M15LS-rot. The left plot shows the fluctuating total momentum in the direction of the 
polar (rotation) axis for all gas on the grid, the right plot displays the corresponding displacement of the center of mass from the grid center 
as function of post-bounce time. Until near the end of the computed evolution, this displacement is far less than the radial width of the 
innermost grid zone (0.3 km). Only when the postshock gas gains momentum as the explosion takes off, which happens with more strength 
in the southern hemisphere and thus corresponds to a growing negative momentum value, the z-displacement exhibits a clear trend and 
becomes slightly larger than the innermost grid zone. 
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Fig. 20. — Left: Test simulations with artificially excited core g-mode oscillations of dipole {I = 1) character in the neutron star 
(here in the case of Model M15HW-2D). Two different amplitudes were considered for an initially imposed velocity field, 5 X lO'^cm/s 
and 2 X lO^cm/s, corresponding to a factor of 16 different kinetic energies (as indicated in the plot). The time evolution of different 
spherical harmonics components of {vr{r,8) — {vr{r, 6)) g) is shown at a radius of 10km. The instigated oscillation remains basically of 
dipole character. The clear presence of many cycles of the oscillation demonstrates the ability of our numerical code to follow such gravity 
waves, if they are instigated. In the lower panel three overlapping curves are displayed. The dashed (hardly visible) black line gives the 
result of a 2D simulation with our standard central ID-core of <1.7km radius, the black solid line is the result with a ID-core of ~0.8km 
radius, and the red bold line shows a simulation in which the whole star down to the center was computed in 2D. Neither the frequency nor 
the amplitude nor the damping behavior are affected by the spherically symmetric treatment of a central region when this ID core is as 
small as chosen. Slightly faster damping and a slowly evolving frequency shift are observed when the ID-core is increased to 3.0 km radius 
(dashed cyan line). Right: The amplitude of the 1 = 1 mode shown in the lower left panel (with a ID-core radius of 0.8 km) is displayed 
as a function of time and radius. The upper panel displays the fractional pressure variations (in percent), {P{r,9) — (Pi^jd)) g)/ {P)g, the 
lower panel the velocity variations normalized by 10^ cm, {vr{r, 6) — {vr{r, 6))g)/i-0^ cm. Note that in both panels the range of values on 
the color bar is limited (cutting off the true maxima and minima of the amplitudes) for visualizing activity in all regions. Interior to about 
10 km the core oscillates with twice the frequency of the mantle outside of r fa 25 km. In the intermediate, convective layer the gravity 
waves are damped and a frequency change happens. At radii r > 50 km one can see the presence of gravity waves produced by the violent 
and chaotic influence of convective overturn and SASI activity in the layers beween the neutron star and the shock. 



